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ABSTRACT 

This  thesis  presents  several  computer  models  used  in  the  design  and  analysis  of  a 
Hit-to-Kill  Rocket  Guidance  Kit  (HRGK).  The  HRGK — ^proposed  as  an  inexpensive  add¬ 
on  kit — ^has  the  potential  of  converting  unguided  2.75”  diameter  rockets  into  precision 
weapons  against  non-tank  targets.  A  Naval  Postgraduate  School  design  team  recently 
participated  in  a  nation-wide  graduate  student  competition  for  the  design  of  such  a  kit. 
The  design  and  analysis  process  led  the  author  to  develop  and  use  various  computer 
models  and  simulations.  This  thesis  documents  three  distinct  types  of  computer  models 
found  useful  in  the  design. 

The  first,  operational  effectiveness  modeling,  established  the  cost  effectiveness  of 
the  NPS  HRGK.  The  second  was  related  to  the  preliminary  sizing  of  various  design 
aspects — ensuring  the  proper  flow-down  of  system  requirements  into  design 
specifications.  The  third  was  a  six-degree  of  freedom  (6DOF)  simulation,  developed  to 
perform  detailed  analyses  of  the  HRGK’s  performance. 

Although  the  models  presented  in  this  thesis  pertain  to  the  HRGK,  the  basic 
principles  apply  to  the  design  or  evaluation  of  other  missile  systems,  and  this  thesis 
provides  general  insights  regarding  the  benefits  and  limitations  of  computer  modeling  in 
missile  design. 


V 


r 


VI 


TABLE  OF  CONTENTS 


l.  INTRODUCTION . 1 

A.  BACKGROUND . 1 

B.  NFS  HRGK  DESIGN . 2 

C.  THESIS  SCOPE  AND  OUTLINE . 4 

n.  COST  OF  OPERATIONAL  EFFECTIVENESS  ANALYSES . 7 

A.  SIMPLIFYING  ASSUMPTIONS . 7 

B.  COST  PER  KILL . .10 

C.  KILLS  PER  SORTIE . 15 

D.  COST  PER  MISSION . 15 

E.  OBSERVATIONS . 20 

m.  PRELIMINARY  SIZING  ANALYSES . 23 

A.  LENGTH,  WEIGHT,  AND  CENTER  OF  GRAVITY . 23 

B.  MANEUVERABILITY . 26 

C.  SEEKER  FIELD  OF  VIEW . 39 

D.  COMMENTS  AND  OBSERVATIONS . 46 

IV.  SIX-DEGREE  OF  FREEDOM  SIMULATION . 51 

A.  SIMULATION  OVERVIEW . ! . 51 

B.  SIMULATION  COMPONENTS . 54 

C.  EQUATIONS  OF  MOTION . 55 

D.  NPS  HRGK  RESULTS . 58 

E.  OBSERVATIONS . 60 

V.  CONCLUSIONS  AND  RECOMMENDATIONS . 61 

APPENDIX  A.  MISSILE  AERODYNAMICS  NOMENCLATURE . . : . 63 

APPENDIX  B.  NOTES  ON  MATLAB  AND  SIMULINK . 65 

APPENDIX  C.  MATLAB  CODE  LISTINGS . 67 

LIST  OF  REFERENCES . . 89 

INITIAL  DISTRIBUTION  LIST . 91 


LIST  OF  FIGURES 


1.  Sample  Normal  Plane  View  for  a  12  mil,  One-Dimensional  Dispersion  at  5  km 

Downrange . 3 

2.  Schematic  View  of  the  NPS  HRGK  Design . 5 

3 .  Projected  Area  for  Light  Armored  Vehicles . 9 

4.  Single- Shot  and  Number  of  Rockets  Required  per  Target  as  a  Function  of 

CEP . 13 

5 .  Relative  Cost  per  Kill  against  Various  Targets . 14 

6 .  AH-64  Apache  Attack  Helicopter  Sorties  Required  to  Achieve  Mission  Objective.  . .  19 

7 .  Number  of  Weapons  Fired  to  Achieve  Mission  Objective . 19 

8 .  Number  of  Targets  Killed  by  Category,  Weapon,  and  Loadout  Option . 20 

9.  Physical  Size  and  Center  of  Gravity  Design  Space . 25 

1 0.  Graphical  Representation  of  Maneuverability  Model  Inputs . 27 

1 1 .  Closed-Form  Solution  for  Maneuverability  Specifications  with  a  Stationary  Target. .  29 

12.  Moving  Target  Scenario  for  Maneuverability  Analysis . 31 

13.  Maximum  Turn  Rate  Required  to  Hit  Target  (95%  Probability) . 33 

14.  Minimum  Turn  Radius  Required  to  Hit  Target  (95%  Probability) . 33 

15.  Plane- View,  Free-Body  Diagram  for  Steady-Tum  Maneuver. . 35 

16.  Required  Thrust  and  Angle  of  Attack  for  Maximum  Turn  Rate  or  Minimum  Turn 

Radius . 38 

17.  Required  Field  of  View.  Results  for  95%  Probability  of  Target  Staying  in  Field  of 

View  with  Worst  Case  Target  Motion . 40 

1 8 .  Target  Distribution  Model . 42 

19.  Illustrative  Case  for  Probabilistic  Seeker  Field  of  View  Analysis . 43 

20.  Possible  Target  Locations . 44 

2 1 .  Alternative  Method  for  Computing  Target  Distribution  Area  within  Seeker  Field  of 

View  Limits . 45 

22.  Probability  Density  Functions,  (a)  Rocket  cross-track  error,  (b)  Target  cross-track 

position,  and  (c)  Probability  of  acquisition . 45 

23.  Probability  of  Target  in  Field  of  View  for  Maximum  Launch  Range  (6  km) . 47 

24.  Sensitivity  Analysis  of  the  Steady-Tum  Maneuver  Model  for  Thmster  Sizing . 49 

25.  HRGK  6DOF  Simulation  Overall  Block  Diagram . 52 

26.  Detailed  View  of  the  Rocket  Dynamics  Block . 53 

27 .  Detailed  View  of  the  Parameters  and  Coefficients  Block . 53 

2  8 .  Future  6DOF  Architecture . 54 

29.  HRGK-Equipped  Rocket  Ballistic  Trajectories . 59 

30.  Sustained  and  Pulsed  Thruster  Maneuvers . 59 


IX 


r 


X 


LIST  OF  TABLES 


1.  Measures  of  Merit  and  Their  Associated  Objectives  and  Requirements . 3 

2.  COEA  Target  Set  with  Assumed  Dimensions,  Desired  and  . 9 

3.  Single-Shot  Probability  of  Kill  Computer  Model  Inputs . 11 

4.  Cost  per  Kill . 14 

5.  Kills  per  AH-64  Apache  Sortie  against  Light  Armor  Targets . 15 

6.  Weight,  Length,  and  Center  of  Gravity  Model  Inputs . 24 

7.  Maneuverability  Model  Inputs . 26 

8.  Steady-Tum  Maneuver  Computer  Model  Inputs . 36 

9.  Seeker  Field  of  View  Analysis  Inputs . 41 


r 


xii 


ACKNOWLEDGMENT 


I  am  grateful  for  the  support  of  my  thesis  advisor,  Dr.  Conrad  Newberry,  and  the 
second  reader,  Dr.  Lou  Schmidt,  for  their  help  and  guidance.  Dr.  Robert  Ball,  Aeronautics 
and  Astronautics  Department,  and  Major  Vince  Tobin  and  Captain  Pat  Mason,  U.S.  Army, 
were  also  very  helpful  with  the  cost  of  operational  effectiveness  portion  of  this  thesis.  I  am 
very  appreciative  of  the  sacrifices  of  Major  Boaz  Pomerantz,  Israeli  Air  Force;  Major 
Silvino  L.  Silva,  Brazilian  Air  Force;  and  Lieutenant  Nigel  A.  Nurse,  U.S.  Navy  who 
served  on  the  NPS  Hit-to-Kill  Rocket  Guidance  Kit  Design  Team. 

I  am  indeed  indebted  to  my  management  at  the  Naval  Air  Warfare  Center  Weapons 
Division  and  the  NAWCWPNS  Long-Term  Fellowship  Committee  for  their  continued 
support  of  my  studies.  I  wish  to  specifically  acknowledge  Tomma  Bersie,  Doug  Savage, 
John  Freeman,  Lee  Gilbert,  Paul  Homer,  and  Ron  Derr  for  allowing  me  this  time  at  the 
Naval  Postgraduate  School.  Special  thanks  also  go  to  Ted  Fincher  and  Ken  Morton  for 
inspiring  me  to  pursue  a  graduate  education. 

But,  my  most  tender  and  heart-felt  gratitude  is  reserved  for  my  sweet  wife  Diane 
and  our  dear  children:  Andrea,  Jared,  Nathan,  Melissa,  Marsha,  Jennifer,  and  Kayla. 
Without  their  love  and  sacrifice,  this  thesis  would  not  have  been  possible. 


xiii 


r 


XIV 


I.  INTRODUCTION 


During  the  first  half  of  the  1996-1997  Naval  Postgraduate  School  (NPS)  academic 
year,  a  team  of  aeronautical  engineering  graduate  students  was  assembled  to  perform  the 
eonceptual  design  of  a  hit-to-kill  guided  rocket  kit  (HRGK).  (The  kit  attaches  to  existing 
unguided  rockets  to  give  them  precision  strike  capability.)  As  part  of  the  design  team,  the 
author  developed  several  computer  codes  or  models  to  facilitate  the  design  and  analysis  of 
the  HRGK.  Since  the  completion  of  the  design  project,  many  of  these  codes  have  been 
expanded  or  refined. 

This  thesis  documents  the  theory  used  in  the  development  of  the  improved 
computer  models.  It  also  describes  the  process  of  using  the  models  and  presents  the 
modeling  results  for  the  NPS  HRGK  design.  These  descriptions  provide  an  example  of  the 
type  of  modeling  that  can  be  used  in  the  design  and  analysis  of  missiles. 

This  introductory  chapter  provides  background  information  about  the  HRGK  and 
pertinent  details  of  the  NPS  HRGK  design.  It  also  outlines  the  remaining  chapters  of  the 
thesis. 


A.  BACKGROUND 

In  the  mid  1990s,  US  Army  Aviation  identified  the  need  for  a  low  cost,  precision 
kill  weapon  system  for  use  against  soft  or  light-armor  targets  (such  as  trucks,  armored 
personnel  carriers,  artillery,  air  defense  systems,  command  posts,  amphibious  landing 
vehicles,  or  patrol  boats).  The  objective  of  such  a  system  would  be  to  minimize  the  cost  per 
kill  for  these  types  of  targets.  The  2.75-inch  rocket  (Mark-66)  with  a  unitary,  high- 
explosive  (HE)  warhead  (M151)  and  a  point  detonation  fiize  (M423)  was  subsequently 
identified  as  a  promising  candidate  for  this  low-cost  mission. 

These  unguided  rockets  cost  only  about  $1,000  per  unit,  and  large  inventories  are 
already  available  to  the  military  services.  The  2.75-inch  rockets  are  launched  from 
helicopters,  fixed-wing  aircraft,  surface  combat  vehicles,  and  naval  vessels.  U.S.  Army 
Helicopters  carry  the  rockets  in  either  a  7-  or  19-tube  launcher  pod  as  part  of  the  Hydra-70 
weapon  system.  The  Hydra-70  system  includes  the  weapon  (rocket,  fuze,  and  warhead), 
launch  pod,  and  an  armament  management  system.  With  this  system,  the  rockets  can  be 
fired  in  single-shot  mode  or  in  various  sized  ripple  or  salvo  shots.  The  rockets  can  carry 
various  warheads  or  payloads  (including  HE,  kinetic  energy  flechettes,  multi-purpose 
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submunitions,  illumination  flare,  marking  smoke,  and  chaff)  with  various  fuzing  options 
(including  remotely  set  timed  burst  and  super-quick  or  delayed  contact  fuzes). 

The  Mark-66  rocket  with  warhead  and  fuze  is  spin  stabilized  to  provide  some 
degree  of  accuracy.  The  rocket  uses  a  fluted  single  throat  nozzle  to  achieve  approximately  a 
10  Hz  spin  by  the  time  it  exits  the  launcher.  By  motor  burnout,  the  spin  rate  has  increased 
to  approximately  35  Hz.  After  motor  bum  out,  the  rocket’s  wrap-around  fins  de-spin  the 
projectile  and  cause  a  significant  spin  rate  (15-20  Hz)  in  the  reverse  direction.  The  spinning 
characteristics  of  the  rocket  played  a  significant  role  in  the  NFS  design  as  will  be  discussed 
in  the  following  section. 

The  baseline  rocket  configuration  chosen  for  the  kit  (HE  warhead  and  super-quick 
contact  fiize)  is  an  area  suppression  weapon.  The  unguided  rocket  has  a  total  12 
milliradians  (mill),  one-dimensional  dispersion,  meaning  that  for  every  thousand  meters  of 
down  range  flyout  a  12  m  error  results  (1  a  or  68%).  For  example,  a  nominal  5  km  range 
launch  would  result  in  a  normal-plane’,  radial  miss  of  up  to  60  m  (5  x  12) — with  39% 
probability  (see  Figure  1).  This  level  of  accuracy  is  far  from  that  required  for  single-shot, 
precision  kills.  To  achieve  point  target  accuracy,  the  rocket  needs  to  be  fitted  with  some 
type  of  guidance  and  control  kit.  The  conceptual  NFS  HRGK,  described  in  the  next 
section,  was  designed  to  serve  this  purpose. 

B.  NFS  HRGK  DESIGN 

The  NFS  HRGK  was  designed  to  meet  certain  requirements.  The  following 
subsections  summarize  these  requirements  and  the  pertinent  aspects  of  the  NFS  HRGK 
design. 

1.  HRGK  Design  SpeciHcations 

The  NFS  HRGK  design  requirements  were  prepared  by  the  American  Institute  of 
Aeronautics  and  Astronautics  (AIAA)  Missile  Systems  Technical  Committee  and  were  used 
for  the  1996/1997  Graduate  Team  Missile  Design  Competition  sponsored  by  the  AIAA  and 
Northrop  Grumman  Corporation  [Ref.  1].  Table  1  summarizes  the  HRGK’s  measures  of 
merit  with  their  associated  objectives  and  pertinent  requirements. 


’  For  the  purposes  of  this  thesis,  all  accuracies  are  in  terms  of  the  plane  normal  to  the  rocket  trajectory 
rather  than  in  the  ground  plane.  This  convention  is  typical  for  guided  missiles. 
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Figure  1.  Sample  Normal  Plane  View  for  a  12  mil,  One-Dimensional 
Dispersion  at  5  km  Downrange.  [(12  m  /  1000  m)5000  m  =  60  m  (la)]. 


Table  1 .  Measures  of  Merit  and  Their  Associated  Objectives  and  Requirements. 


Measure  of  Merit  and  Objective 

Specific  Requirements 

Cost 

Maximize  Units  Bought 

•  $10,{X)0  Avg.  Unit  Cost  (5,000  Unit  LRIP") 

•  Maximum  $1  Billion  Total  5  Year  Production  Buy 

Current  Systems  Computability 
Maximize  Operability  & 

Reduce  Life  Cycle  Cost 

•  Existing  Rocket,  Fuze,  Warhead,  &  Launch  Pod 

—  72  inch  total  system  maximum  length 

—  loaded  pod  center  of  gravity  between  lugs 

•  Laser  Designator  Compatibility 

Accuracy 

Provide  Hit-to-Kill  Accuracy 

•  0.5  m  CEP 

•  60  mph  target  speed 

Weight 

Maximize  Kits  Deliverable  via  Airlift 

•  30  pound  Total  System  Maximum  Weight 

Maximum  Range 

Maxinndze  Platform  Standoff  Survivability 

•  6  km  (or  More)  Maximum  Range 
(Maneuverability  at  Maximum  Range) 

Minimum  Range 

Maximize  Launch  Acceptable  Region 

•  1  km  (or  Less)  Minimum  Range 

(Guidance,  Fuzing,  and  Launch  Transients) 

Adverse  Weather  Performance 

Maximize  Operability  and  Combat  Flexibility 

•  99%  Worldwide  Weather  Capable 

Time  to  Target 

Minimize  Time  (Improve  Survivability) 

a  T _ i  v-k  1 

•  16  seconds  or  Less  to  5  km  Range 

"  Lx)w  Rate  Initial  Production 


3 


2.  NPS  Design  Highlights 

The  following  paragraphs  briefly  describe  the  major  features  of  the  NPS  HRGK 
design.  The  design  team’s  proposal  gives  a  more  detailed  description  of  the  design  [Ref. 
2]. 

The  HRGK  attaches  over  the  fuze  at  the  front  end  of  a  2.75  inch  rocket.  It  utilizes  a 
strapped-down,  four-quadrant  silicon  detector  for  semi-active  laser  spot  homing.  The 
HRGK  employs  three  hot  gas  reaction  jets  (thrusters)  to  control  the  rocket’s  attitude  and 
trajectory.  A  solid  propellant  burning  in  an  insulated  pressure  chamber  generates  the 
thruster  gas.  The  kit  uses  a  look-up  table  scheme  to  control  the  thrusters  based  solely  on  the 
seeker  outputs.  Fast-cycle  solenoid  valves  open  and  close  the  “on-off’  thrusters.  The 
thruster  nozzles  are  inter-spaced  between  three  fixed,  wrap-around  canards.  The  canards 
provide  increased  maneuverability  (by  reducing  the  rocket’s  static  stability)  and  help 
constrain  any  jet  and  air  stream  interaction  effects.  The  canards  also  contribute  significantly 
to  the  HRGK’s  roll  stability. 

Thrust  and  radial  needle  rollers  effectively  isolate  the  NPS  HRGK  from  the 
spinning  motion  of  the  rocket.  Simulations  show  the  kit’s  roll  rate  to  be  less  than  9  degrees 
per  second  after  the  first  seconds  of  flight.  The  kit’s  low  spin  rate  and  the  simplistic  pursuit 
navigation  scheme  eliminate  the  need  for  inertial  sensors  and  greatly  reduce  the  kit’s  signal 
processing  requirements. 

As  shown  in  Figure  2,  the  basic  kit  is  19.0  inches  long.  It  has  a  diameter  of  2.75 
inches  (with  canards  wrapped  in)  and  weighs  7.0  pounds.  The  NPS  HRGK  attaches  to  the 
2.75  inch  rocket  between  the  Ml 51  high-explosive  warhead  and  M423  point  detonation 
fuze.  The  all-up  round  is  68.75  inches  long  and  weighs  30  lb.  The  guided  rocket  is 
compatible  with  the  Hydra-70  rocket  system. 

C.  THESIS  SCOPE  AND  OUTLINE 

The  NPS  HRGK  design  was  achieved  and  verified  with  the  aid  of  several  original 
computer  codes  and  models.  The  codes  and  models  developed  by  the  author  and  described 
in  this  thesis  fall  into  the  following  three  major  categories: 

•  cost  of  operational  effectiveness  analyses, 

•  preliminary  sizing  analyses,  and 

•  six-degree  of  freedom  simulation. 
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Figure  2.  Schematic  View  of  the  NFS  HRGK  Design.  Top  drawing  shows  HRGK  in  place  on  Mark-66  Rocket  with  warhead  and 
fuze.  The  in-flight  ranges  of  the  guided  rocket’s  center  of  gravity  and  center  of  pressure  are  also  indicated  on  the  top  drawing. 


The  following  three  chapters  provide  descriptions  of  the  modeling  processes  in  the 
respective  categories.  Each  of  these  chapters  (1)  includes  a  description  of  the  modehng 
methodology,  (2)  gives  example  results  specific  to  the  NFS  HRGK  design,  and  (3) 
discusses  observations  regarding  the  general  use  of  the  models.  The  final  thesis  chapter 
summarizes  these  conclusions  and  provides  recommendations  for  additional  work  with  the 
computer  models.  Appendices  that  contain  the  computer  code  listings  and  additional 
background  have  also  been  included.  The  appendices  include  a  brief  description  of  missile 
aerodynamics  nomenclature  (Appendix  A)  and  a  discussion  concerning  the  use  of  the 
computer  software  Matlab®  and  Simulink®  (Appendix  B). 

In  general,  this  thesis  is  not  intended  to  be  read  from  cover  to  cover.  Instead  this 
thesis  is  offered  as  a  reference  to  possible  analysis  types  and  modeling  methods  for  use  in 
future  conceptual  design  of  missile  systems. 


®  Matlab  and  Simulink  are  registered  trademarks  of  The  MathWorks,  Inc. 
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II.  COST  OF  OPERATIONAL  EFFECTIVENESS  ANALYSES 


Cost  of  operational  effectiveness  analysis  (COEA)  must  be  an  important  part  of  the 
design  process.  First,  COE  As  provide  justification  for  the  design  based  on  operational 
effectiveness  cost  savings.  But,  just  as  importantly,  COEA-type  analyses  provide  the 
designer  with  focus  and  insight  into  the  truly  important  aspects  of  the  design.  Basic  COEA 
models  with  many  simplifying  assumptions  can  be  used  early  in  the  process  to  focus 
design  decisions.  As  the  design  matures,  refined  models  can  be  used  to  evaluate  major 
design  iterations. 

The  COEA  modeling  developed  for  the  NFS  HRGK  provided  the  justification  for  a 
HRGK-lUce  system.  The  overall  objective  of  the  HRGK  was  to  minimize  the  total  cost  per 
kill  against  appropriate  targets.  From  the  analysis,  guided  rockets  showed  a  decreased  cost 
per  target  kill  and  increased  target  kills  per  sortie  when  compared  to  alternative  weapon 
systems.  Additional  modeling  showed  that  guided  rockets  resulted  in  a  reduction  of  sorties 
and  weapon  costs  required  to  achieve  a  specific  mission  objective.  The  final  NFS  HRGK 
was  an  attempt  to  provide  the  lowest  possible  cost  per  kill  design. 

The  COEA  modeling  discussed  in  this  chapter  provides  several  different  levels  of 
results.  The  low  level  results  of  the  modeling  include  single-shot  probabilities  of  kill  and 
the  number  of  weapons  required  per  kill  as  a  function  of  circular  error  probable  (CEF). 
These  values  are  then  used  to  derive  the  cost  per  kill  and  the  sorties  per  mission  results. 
The  number  of  weapons  required  per  kill  results  also  provide  some  additional  insights  that 
was  critical  to  the  NFS  HRGK  design. 

The  following  section  describes  the  simplifying  assumptions  used  in  the  NFS 
HRGK  COEA  models.  Subsequent  sections  provide  more  detailed  descriptions  of  the 
model  methodologies  and  the  results.  The  chapter  ends  with  a  discussion  of  some 
observations  concerning  the  use  of  COEA  modeling  in  design.  (Listings  of  the  computer 
codes  used  in  the  NFS  HRGK  COEA  are  included  in  Appendix  C.) 

A.  SIMPLIFYING  ASSUMPTIONS 

Several  simplifying  assumptions  were  made  to  facilitate  the  early  NFS  HRGK 
COEA  modeling.  These  assumptions  included  defining  guided  and  unguided  weapon 
trajectories  and  costs,  establishing  a  target  set  with  desired  probabilities  of  kill  (P^)  for  each 
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target,  and  estimating  each  target’s  vulnerability  in  terms  of  the  conditional  probability  of 
kill  given  a  hit  {Pan)- 

1.  Weapon  Accuracies 

The  accuracies  for  the  unguided  rocket  were  taken  as  a  function  of  launch  range  and 
were  based  on  a  12  milliradian  accuracy  [Ref.  1].  As  a  baseline  for  the  comparison,  a  mid¬ 
range,  3-km  shot  was  used  with  a  resulting  42.4  m  normal-plane  CEP  (36  m  1- 
dimensional,  standard  deviation  multiplied  by  1.1774 — see  Figure  1,  for  example).  For  the 
guided  rockets,  three  different  normal-plane  CEP  accuracies  was  assumed  for  comparison 
purposes — 0.5  m,  1.0  m,  1.5  m,  and  3.0  m. 

For  both  the  unguided  and  guided  rockets,  a  20°  terminal  dive  angle  (relative  to  the 
horizontal)  was  assumed  as  a  best  estimate'.  The  shallow  dive  caused  an  increase  in  the 
rocket’s  ground  plane  down-track  error  but  did  not  affect  the  cross-track  error.  These 
effects  were  accounted  for  implicitly  in  determining  the  probability  of  hitting  the  target. 

For  simple  comparison  purposes,  the  AGM-114  laser  Hellfire’s  accuracy  was 
assumed  to  be  such  that  if  it  was  shot  at  a  target  it  would  achieve  the  desired  killed. 

2.  Target  Set 

The  effectiveness  analysis  used  four  representative  targets  (mobile  air  defense  unit, 
armored  personnel  carrier,  support  vehicle,  and  patrol  boat).  These  targets  were  each 
assumed  to  be  simple  blocks  of  appropriate  length,  width,  and  height  dimensions  (see 
Figure  3).  Table  2  lists  the  four  targets  with  their  associated  dimensions  [Ref.  3,  4,  5,  and 
6].  The  dimensions  are  representative  for  targets  found  in  the  references.  Table  2  also  lists 
a  desired  Pfr  for  each  of  the  targets.  These  desired  probabilities  were  assumed,  based  on  the 
typical  importance  of  the  target. 

The  probability  of  each  target  being  hit  (P^)  was  computed  as  the  joint  probability 
of  the  rocket  being  within  both  the  top-to-bottom  and  the  side-to-side  dimensions  of  the 
target.  The  projected  dimensions  of  the  targets  (normal  to  the  flight  trajectory)  were 
computed  for  both  a  head-on  and  a  broad-side  attack.  The  two  attack  approaches  gave  the 
extremes  in  the  target’s  presented  area  as  shown  in  Figure  3.  The  P^’s  for  the  two  extreme 
cases  were  computed  and  then  equal-weight  averaged  as  an  estimate  of  the  for  uniformly 
distributed  attack  approaches  (0°-360°). 


*  Later  analysis  showed  that  unguided  rockets  may  have  even  more  shallow  dive  angles,  making  them 
even  more  inaccurate  in  the  ground-plane. 
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Figure  3.  Projected  Area  for  Light  Armored  Vehicles  (head-on, 
45°,  and  broad-side  attacks  with  20°  dive  angle). 


Table  2.  COEA  Target  Set  with  Assumed  Dimensions,  Desired  P,^,  and  P^. 


Target  Description 

Dimensions 

L  X  W  X  H 

Desired 

Mobile  Air  Defense  Unit  (Light  Armor) 

5  X  2.5  X  2  m 

95% 

90% 

Armored  Personnel  Carrier  (APC) 

5  X  2.5  X  2  m 

60% 

80% 

Support  Vehicle  (Large  Truck) 

8  X  2.5  X  2.5  m 

50% 

100% 

Patrol  Boat 

14  X  4.5  X  1.5  m 

75% 

70% 

3.  Target  Vulnerabilities 

For  each  target  a  “cookie-cutter”  model  [Ref.  7]  was  used  in  determining  the  single¬ 
shot  Pf.  With  the  cookie  cutter  model,  the  target  suffers  a  kill  (according  to  an 

assumed  if  and  only  if  it  is  actually  hit.  “Near  misses”  do  not  provide  any  P^.  Thus 

the  value  for  Pf,  was  the  simple  product  of  P„^  and  P^  .  The  assumed  values  of  Pj^,^  are 
best  estimates  and  are  shown  in  the  last  column  of  Table  2.  These  high  values  of  Pjf,^ 
reflect  the  facts  that  (1)  the  targets  were  relatively  soft,  (2)  the  rocket  is  actually  impacting 
the  target  to  score  a  hit,  and  (3)  the  desired  type  of  kill  (firepower,  mobility,  catastrophic, 
etc.),  which  was  not  specified,  may  only  require  minor  damage  to  the  target. 

4.  Weapon  Costs  and  Additional  Pj^’s 

The  per  unit  cost  of  the  guided  rocket  concepts  was  assumed  to  be  $10,000 — ^the 
upper  threshold  cost  from  the  AIAA  design  contest,  and  the  unguided,  unitary  warhead 
round  (M151)  was  assumed  to  cost  $1,000  [Ref.  1].  A  multipurpose  submunitions 
(MPSM)  round  was  assumed  to  cost  $2,000  and  have  a  10%  Pj^  against  non-tank  targets. 
The  Hellfire  missile  was  assumed  to  have  a  per  unit  cost  of  $40,0(X)  [Ref.  8]  and  a  single¬ 
shot  kill  capability  against  all  four  representative  targets  (P^  =  80%  for  tanks  and  P^  = 
100%  for  all  others). 
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5.  Additional  Assumptions 

Additionally,  the  assumption  was  made  that  all  weapons  functioned  during  launch, 
flyout,  and  fuzing;  with  100%  reliability.  The  aimpoint  was  assumed  to  be  the  centroid  of 
the  target  block  without  any  bias  errors  in  designation  or  flyout.  It  was  also  assumed  that 
all  rocket  accuracy  errors  had  a  bivariate  normal  distribution  with  independent,  equi- 
variance  (a^)  distributions  along  any  two  orthogonal  axes  normal  to  the  flyout. 
Furthermore,  the  errors  were  assumed  independent  from  rocket  to  rocket.  Thus,  the  CEP 
was  1.1774  a  and  the  Pj^  for  n  multiple  rockets  was  1  -  i}-Pic single shoT- 

The  simplifying  assumptions  discussed  in  this  section  are  reasonably  valid  under 
many  conditions.  The  last  section  of  this  chapter  (lI.E.  Observations)  includes  a  brief 
discussion  on  the  affect  of  some  of  the  assumptions. 

B.  COST  PER  KILL 

The  overall  objective  of  the  HRGK  was  to  minimize  the  total  cost  per  kill  against 
appropriate  targets.  To  determine  the  cost  per  kill,  the  probability  of  kill  per  single  weapon 
had  to  be  established  using  the  assumptions  discussed  above.  This  single-shot  probability 
of  kill,  with  the  additional  assumptions  above,  could  then  be  used  to  determine  the 
number  of  weapons  required  to  achieve  the  desired  Pf..  From  the  number  of  weapons 
required  per  kill  and  the  assumed  cost  per  weapon,  the  estimated  cost  per  kill  could  then  be 
established.  The  process  for  determining  the  cost  per  kill  and  the  results  for  the  HRGK 
analysis  are  discussed  in  the  following  subsections. 

1.  Methodology 

The  methodology  for  determining  the  cost  per  kill  involves  determining  the  P^^  for 
each  weapon  against  each  target  and  then  determining  the  number  of  weapons  required  to 
achieve  the  desired  P^.  A  computer  program  was  written  in  Matlab  to  compute  the  P^^^ 
and  number  of  weapons  required  per  kill  as  a  function  of  CEP.  This  program,  which  also 
computes  a  relative  cost  per  kill  tradeoff  plot,  is  described  in  the  following  paragraphs. 

a.  Model  Inputs 

Only  a  small  number  of  inputs  are  required  for  the  single-shot  probability  of 
kill  computer  model.  These  inputs  can  be  divided  into  descriptions  of  the  target  and  of  the 
weapon  trajectory.  Table  3  lists  these  two  types  of  required  inputs. 
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Table  3.  Single-Shot  Probability  of  Kill  Computer  Model  Inputs. 


Input  Name 

Input  Description 

Target 

L,  H,  and  W 

Effective  Target  Length,  Height,  and  Width  (unit  of  length) 

Pkh 

Probability  of  Kill  given  a  Hit  within  the  Effective  Target  Dimensions 

PkD 

Desired  Probability  of  Kill  for  the  Target 

Weapon  Trajectory 

dive 

Dive  Angle  Relative  to  Horizontal  (degrees) — 9 

CEP_min 

Minimum  CEP  to  Be  Used  in  Computations  (unit  of  length) 

CEP__max 

Maximum  CEP  to  Be  Used  in  Computations  (unit  of  length) 

b.  Single-Shot  Probability  of  Kill 

The  first  step  in  estimating  the  single-shot  probability  of  kill  is  to  estabhsh 
the  end-on  and  broad-side  projected  edges  of  the  target  normal  to  the  weapon  trajectory. 
The  side-to-side  dimensions  for  the  two  cases  are  simply  Xj  =W  and  X2  =  L,  for  the  end- 
on  and  broad-side  cases,  respectively.  The  top-to-bottom  dimensions  are  given  by  the 
following  equations: 

=  Hcos  0  Lsin  0  (end-on) 

Y2  =  H  cos  d  -^W  sin  d  (broad-side) 

Next,  the  probability  of  hit,  P„,  was  determined  by  estimating  the 

probability  of  the  weapon  striking  within  the  projected  edges  of  the  target — ^the  joint 

X.  Y- 

probability  of  being  within  ±— i-  and  ±—  of  the  aimpoint,  where  i  =  1  or,  2  depending  on 

the  approach  case.  Given  the  assumed,  independent  bivariate  nature  of  the  weapon  errors, 
the  joint  probability  could  be  computed  as  the  product  of  the  two — cross-track  (Z- 
direction)  and  elevation  (T-direction) — ^probabilities.  The  erf  function  in  Matlab  is 
defined  as  the  following: 


so  the  joint  probability  described  above  is  given  by  the  following  Matlab  expression 
where  the  argument  to  erf  is  normalized  by  : 
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y  y  ^ 


=  erf 


X  1 


erf 


2  -Jltr)  U^/2<T 


and  where  /  =  1  or  2  and  cris  the  CEP/1.1114.  The  total  P„  is  then  the  equal  weighted 
average  of  P„^  and  Pfj^.  In  other  words,  P„  =  {P„  +  P„  )/2. 

With  the  probability  of  hit,  P^,  and  the  input  Pg,„,  the  single-shot 
probability  of  kill  is  estimated  as  the  product  of  the  two. 

P  =  P  P 

^Kss  ^KIH 


c.  Number  of  Weapons  per  Kill 

With  the  Pf^^  and  the  assumption  that  flyout  errors  for  each  weapon  shot  are 
unbiased  and  independent,  the  number  of  weapons  required  to  achieve  a  desired  Pg^  can  be 
calculated. 

The  probability  of  the  target  surviving  a  single  shot  is  given  as  the 
complement  of  l-P^ss-  Th®  probability  of  surviving  n,  independent  shots  is  that 
quantity  raised  to  the  nth  power,  (l-P^^^^)".  Thus,  the  probability  of  kill  for  n,  independent 
shots  is  the  complement  to  the  latter  probability  of  survival,  namely: 

The  number  of  shots,  n,  to  achieve  a  desired  P^  can  be  determined  by 
solving  the  above  equation  for  n  as  is  shown  in  the  following  equation. 


d.  Accuracy  and  Relative  Cost  Trade-off  Chart 
Knowing  the  number  of  weapons,  n,  required  to  achieve  a  desired 
probability  of  kill  (P^)  against  a  specific  target,  the  relative  weapon  costs  to  give  equal  cost 
per  kill  can  be  computed.  For  example,  if  two  weapons  (A)  with  a  0.5  m  CEP  or  six 
weapons  (5)  with  a  2.0  m  CEP  would  be  required  to  achieve  a  desired  Pj^.,  then  if  the  cost 
of  weapon  B  was  one-third  the  cost  or  less  of  weapon  A,  weapon  B  would  have  a  lower 
cost  per  kill  against  the  specified  target. 

In  the  computer  model,  weapon  costs  are  normalized  relative  to  the  cost  of 
the  minimum  input  CEP,  CEPjnin,  weapon.  Also,  the  number  of  weapons  required  to 
achieve  a  desired  Pj^  are  considered  only  as  integer  values. 
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e.  Cost  per  Kill  Calculations 

The  absolute  cost  per  kill  can  easily  be  computed  from  the  results  of  the 
computer  model  described  above.  The  cost  per  kill  for  a  specified  target  and  weapon  is 
simply  the  integer  number  of  weapons  required  to  achieve  the  desired  against  that  target 
multiplied  by  the  assumed  cost  of  that  weapon. 

2.  Results 

The  results  of  the  computer  model  described  above  are  shown  for  the  targets  and 
weapons  described  above  in  section  II.A.  Simplifying  Assumptions.  The  results  include 
both  the  and  the  cost  per  kill. 

a.  Direct  Model  Outputs 

The  combined  single  shot  probability  of  kill  and  the  number  of  weapons 
required  to  achieve  a  desired  Pg  are  shown  as  a  function  of  CEP  in  Figure  4.  Various 
assumed  accuracies  for  the  unguided  rockets  are  also  marked  on  the  plot  with  their 
corresponding  launch  ranges. 


Figure  4.  Single-Shot  Pg  and  Number  of  Rockets  Required  per  Target  as  a  Function  of  CEP. 
(.Pkss  curves  slope  down  from  upper  left;  number  of  rockets  curves  slope  up  to  upper  right.) 
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Figure  5  shows  the  relative  weapon  costs  for  equal  cost  per  kill  as  a 
function  of  CEP.  In  this  chart  the  weapon  costs  are  relative  to  the  0.5  m  CEP  weapon  cost. 
From  the  chart  it  appears  that  a  weapon  with  a  2  m  CEP  {B)  would  have  an  improved  cost 
per  kill  against  all  targets  if  the  weapon  cost  were  less  than  33%  of  the  0.5  m  CEP  weapon 
(A).  Weapon  C,  at  half  the  cost  of  A,  has  the  highest  cost  per  kill  at  4  m  CEP. 


Figure  5.  Relative  Cost  per  Kill  against  Various  Targets.  Lettered  design  points  (A,  B, 
and  C)  are  shown  for  illustration  only. 


b.  Cost  per  Kill  Results 

From  the  number  of  weapons  required  per  target  kill,  the  costs  per  kill  for 
each  of  the  weapons  can  be  computed.  The  results  for  the  HRGK  analysis  based  on  the 
single  assumed  guided  rocket  cost  are  shown  in  Table  4. 


Table  4.  Cost  per  Kill. 


Target 

Unguided 

Rocket 

Guided  Rocket 

Hellfire 

0.5  m 

1.0  m 

1.5  m 

3.0  m 

Mobile  Air  Defense  Unit 

$2,219k 

$20k 

$20k 

B^fll 

$120k 

$40k 

Light  Armor 

$764k 

$10k 

$l0k 

$50k 

$40k 

Large  Truck 

$294k 

$10k 

$10k 

$10k 

$20k 

$40k 

Patrol  Boat 

$467k 

$20k 

$20k 

$20k 

$40k 

$40k 
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c. 


KILLS  PER  SORTIE 


A  shortcoming  of  the  cost  per  kill  analysis  previously  presented  is  illustrated  in  the 
following  example.  A  weapon  with  a  2.0  m  CEP  (such  as  fi  in  Figure  5)  and  a  quarter  of 
the  cost  of  a  0.5  m  CEP  weapon  has  a  lower  cost  per  kill  than  the  0.5  m  CEP  weapon 
(such  as  A  in  Figure  5).  However,  it  would  take  three  times  as  many  of  the  less  accurate 
(but  inexpensive)  weapons  to  achieve  the  same  number  of  kills  against  mobile  air  defense 
targets.  The  number  of  weapons  available  for  a  given  mission  would  be  limited  by  the 
number  of  launch  platforms  sorties.  The  importance  of  kills  per  sortie  is  completely  ignored 
in  the  cost  per  kill  data.  A  kills  per  sortie  analysis  helps  to  illustrate  the  importance  of 
precision  in  strike  warfare. 

The  kills  per  sortie  analysis  performed  for  the  HRGK  was  based  on  AH-64  Apache 
attack  helicopter  sorties.  The  helicopters  were  assumed  to  have  either  pure  or  mixed 
weapon  payloads  or  loadouts.  The  pure  loadouts  were  either  16  Hellfires  or  76  rockets 
(four  19-tube  rocket  launchers).  The  mixed  loadouts  combined  eight  Hellfires  with  a 
limited  number  of  unitary  warhead  (M151)  rockets  (guided  or  unguided)  and  additional 
MPSM-equipped  rockets.  The  number  of  kills  per  sortie  were  calculated  based  on  the 
data  in  Figure  4  and  the  assumed  P/s  from  the  sub-section  II.A.4.  Weapon  Costs  and 
Additional  PR’s.  The  results  are  tabulated  in  Table  5  for  the  light  armor  target  case. 


Table  5.  Kills  per  AH-64  Apache  Sortie  against  Light  Armor  Targets. 


Weapon 

Maximum  Load 

Mixed  Load 

Hellfire 

16  Hellfires 

16  Kills 

8  Hellfires  &  38  MPSM 

9.90  Kills 

HRGK  (0.5m) 

76  HRGKs 

60.7  Kills 

8  Hellfires,  24  MPSM,  &  14  HRGKs 

20.4  Kills 

HRGK  (1.0m) 

76HRGKS 

52.5  Kills 

8  Hellfires,  24  MPSM,  &  14  HRGKs 

18.9  Kills 

HRGK  (1.5m) 

76  HRGKs 

38.7  Kills 

8  Hellfires,  24  MPSM,  &  14  HRGKs 

16.3  Kills 

HRGK  (3.0m) 

76  HRGKs 

15.1  Kills 

8  Hellfires,  24  MPSM,  &  14  HRGKs 

12.0  Kills 

Unguided 

76  Rockets 

0.09  Kills 

8  Hellfires,  24  MPSM,  &  14  Rockets 

9.22  Kills 

D.  COST  PER  MISSION 

A  final  type  of  modeling  performed  by  the  author  for  the  NFS  HRGK  design 
combines  the  attributes  of  the  cost  per  kill  and  the  kills  per  sortie  analyses.  This  analysis 
assesses  the  cost  (in  terms  of  sorties  and  cost  of  weapons  fired)  to  achieve  a  mission 
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objective.  The  following  subsections  outline  the  analysis  and  present  the  results  for  the 
HRGK. 


1.  Mission  Definition 

The  cost  per  mission  analysis  was  based  on  a  specific  assumed  mission  scenario. 
This  scenario  defined  the  target  matrix  with  a  set  priority  on  the  target  types  as  well  as  a 
mission  objective.  The  scenario  also  defined  the  attacking  helicopter  loadouts.  The  assumed 
scenario  in  defined  in  the  following  paragraphs. 

a.  Target  Matrix 

The  assumed  target  was  an  assembly  area  with  the  following  targets  types 
listed  by  priority  and  with  the  target  quantities  in  parentheses: 

•  Mobile  Air  Defense  Units  (8) 

•  Main  Battle  Tanks  (40) 

•  Light  Armored  Vehicles  (150) 

•  Support  Vehicles  (250) 

The  mission  objective  was  to  destroy  half  of  the  vehicles  in  each  combatant  category. 

b.  Weapon  Loadouts 

The  assumed  initial  strike  force  was  a  flight  of  five  AH-64  Apache  attack 
helicopters.  As  additional  helicopters  were  needed  they  were  added  to  the  attack  force  one 
at  a  time.  Each  Apache  carried  a  nominal  load  of  eight  Hellfires  and  two  19-rocket  launcher 
pods.  Each  pod  (which  could  be  divided  into  zones)  was  loaded  with  twelve  M261  MPSM- 
equipped  rockets  and  seven  guided  or  unguided  M151  unitary  warhead  equipped-rockets. 
[Ref.  9] 

The  limited  number  of  HRGK-equipped  weapons  (14)  per  helicopter  was 
based  on  the  assumption  that  only  six  to  eight  laser-designated  flyouts  could  be  made  from 
a  firing  position  before  the  helicopters  would  be  forced  to  relocate.  With  eight  Hellfires  and 
14  HRGKs,  three  firing  positions  would  be  required  for  the  attack.  This  was  considered  to 
be  the  upper  limit  for  a  deep  strike  mission.  The  MPSM  rockets  would  not  require  laser 
designation,  yet  they  would  still  be  effective  against  clusters  of  light  armor,  material,  or 
personnel.  [Ref,  9] 

2.  Example  Calculation 

The  following  provides  an  example  of  the  process  used  to  model  the  mission 
engagement.  The  1.0  m  CEP  case  is  illustrated. 
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a.  First  Five  Helicopters 

The  weapons  loadout  of  the  first  five  helicopters  consists  of  the  following: 

•  40  Hellfire  Missiles  (each  at  a  tank) 

•  40  X  0.8  =  32  tank  kills 

•  70  HRGK  Rockets  (16  at  air  defense  units — 2  per  unit,  54  at  light  armor) 

•  16 — 8  air  defense  unit  kills  with  =  95%) 

•  54  X  0.6909  “  37.31  light  armor  kills 

•  120  MPSM  Rockets  (1 20  x  0.1  =  12  kills  divided  between  non-tank  targets) 

•  6  light  armor  kills 

•  6  support  vehicle  kills 

In  summary,  after  five  helicopters,  8  air  defense  units  are  kiUed  (exceeding  objective),  32 
tanks  are  killed  (exceeding  objective),  43.31  light  armor  vehicles  are  killed  (31.69  short  of 
the  objective),  and  6  support  vehicles  are  killed  (exceeding  objective). 

b.  Sixth  Helicopter 

The  following  delineates  the  use  of  the  sixth  helicopter’s  weapons: 

•  8  Hellfire  Missiles  (each  at  a  tank) 

•  8  X  0.8  =  6.4  tank  kills 

•  14  HRGK  Rockets  (each  at  light  armor) 

•  14  X  0.6909  =  9.67  light  armor  kills 

•  24  MPSM  Rockets  (24  x  0. 1  =  2.4  kills  divided  between  non-tank  targets) 

•  1.2  light  armor  kills 

•  1.2  support  vehicle  kills 

In  summary,  after  the  sixth  helicopter,  8  air  defense  units  are  killed  (exceeding  objective), 
38.4  tanks  are  killed  (exceeding  objective),  54.18  light  armor  vehicles  are  killed  (20.82 
short  of  the  objective),  and  7.2  support  vehicles  are  killed  (exceeding  objective). 

c.  Seventh  Helicopter 

The  following  is  for  the  seventh  helicopter: 

•  8  Hellfire  Missiles  (2— each  at  a  tank,  6  at  light  armor) 

•  2  X  0.8  =  1.6  tank  kills 

•  6  X  1.0  =  6  light  armor  kills 

•  14  HRGK  Rockets  (each  at  light  armor) 

•  14  X  0.6909  =  9.67  light  armor  kills 

•  24  MPSM  Rockets  (24  x  0. 1  =  2.4  kills  divided  between  non-tank  targets) 

•  1.2  light  armor  kills 

•  1.2  support  vehicle  kills 

After  the  seventh  helicopter,  8  air  defense  units  are  killed  (exceeding  objective),  40  tanks 
are  killed  (exceeding  objective),  71.05  light  armor  vehicles  are  killed  (3.95  short  of  the 
objective),  and  8.4  support  vehicles  are  killed  (exceeding  objective). 
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d.  Subsequent  Helicopters 

For  any  subsequent  helicopters  the  weapons  would  be  used  as  follows: 

•  8  Hellfire  Missiles  (each  at  light  armor) 

•  8  X  1 .0  =  8  light  armor  kills 

•  14  HRGK  Rockets  (each  at  light  armor) 

•  14  X  0.6909  =  9.67  light  armor  kills 

•  24  MPSM  Rockets  (24  x  0.1  =  2.4  kills  divided  between  non-tank  targets) 

•  1.2  light  armor  kills 

•  1.2  support  vehicle  kills 

This  gives  18.87  light  armor  kills  per  sortie;  therefore,  the  number  of  sorties,  after  the 
seventh,  to  reach  the  objective  of  75  light  armor  kills  is  3.95/18.87  =  0.21,  and  the  total 
number  of  sorties  required  for  the  mission  is  7.21. 

e.  Summary 

After  7.21  sorties  the  following  number  of  weapons  have  been  fired  with 
the  indicated  assumed  weapon  cost: 

•  57.67  Hellfires  ($2.3M) 

•  100.9  HRGKs  ($1.0M) 

•  173.0  MPSMs  ($0.35M) 

for  a  total  weapon  cost  of  $3.66M.  The  total  targets  killed  were: 

•  8  air  defense  units 

•  40  tanks 

•  75  light  armor  vehicles 

•  8.65  support  vehicles. 

3.  Results 

Figures  6,  7,  and  8  illustrate  the  results  for  the  HRGK  mission  analysis.  Figure  6 
shows  the  number  of  sorties  required  to  achieve  the  mission  objective.  With  an  assumed 
cost  per  sortie,  these  numbers  could  be  directly  applied  to  the  weapon’s  cost  per  kill.  The 
inclusion  of  a  helicopter  attrition  rate  would  also  make  the  cost  per  kill  more  realistic.  Note 
that  the  1.0  m  CEP  rocket  does  not  require  many  more  sorties  than  does  the  0.5  m  CEP. 

Figure  7  shows  the  number  of  weapons  fired  during  the  mission  for  each  of  the 
different  weapon  cases.  The  total  weapon  cost  in  millions  of  dollars  is  also  indicated  on  the 
graph. 

Figure  8  summarizes  which  weapons  killed  which  targets.  It  is  interesting  to  note 
that  with  accurate  HRGKs,  the  more  expensive  Hellfire  missile  could  be  used  almost 
exclusively  on  tanks.  On  the  other  hand,  when  unguided  rockets  were  used,  nearly  all  the 
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target  kills  were  from  the  Hellfire,  and  the  unguided  unitary  (M151)  rocket  only  registered 
0.242  light  armor  and  0.474  support  vehicle  kills. 


Figure  8.  Number  of  Targets  Killed  by  Category,  Weapon,  and  Loadout  Option. 


E.  OBSERVATIONS 

Simple  COEA  type  modeling  is  essential  to  the  design  process.  Besides  assuring 
that  the  design  can  be  cost  effective,  the  results  from  early  analyses  provide  focus  for  the 
design  team.  COEA  results  can  show  general  trends  that  might  otherwise  go  unnoticed. 
COEA  type  modehng  continues  to  be  useful  throughout  the  design  process  as  refined- 
model  analyses  can  help  in  verifying  design  decisions. 

An  example  of  how  unexpected  information  can  come  from  a  COEA  study  occurred 
in  the  NFS  HRGK  analysis.  Figure  4  (page  13)  shows  a  “knee  in  the  curve”  around  1  m 
CEP,  That  is  the  approximate  CEP  where  the  begin  to  drop  off  significantly  and  the 
number  of  weapons  needed  per  kill  begins  a  steep  increase.  From  the  figure,  it  appears  that 
weapons  achieving  a  CEP  much  smaller  than  1  m  are  not  needed  against  the  assumed 
targets.  Therefore,  as  a  result  of  the  COEA  modeling,  designs  that  achieve  tighter  than  1  m 
CEPs  at  significantly  higher  costs  are  eliminated  from  consideration. 

On  the  other  hand,  COEA  results  can  be  very  sensitive  to  underlying  assumptions. 
The  possibility  of  invalid  assumptions  can  make  COEA  modeling  too  unreliable  for 
confident  design  decisions.  For  example,  small  changes  in  the  low-level  target  set 
assumptions  can  result  in  important  changes  in  values.  These  values,  in  turn,  will 
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impact  sortie  rates.  Likewise,  high-level  assumptions  like  the  mission  definition  used  in  the 
sorties  per  mission  analysis  will  also  impact  the  final  COEA  results.  The  validity  of  the 
assumptions  involved  in  both  low-  and  high-level  models  must  be  continuously  assessed. 

Specific  examples  of  COEA  sensitivities  to  underlying  assumptions  in  the  NFS 
HRGK  models  include  two  opposing  possibilities.  If  the  assumed  target  set  were  expanded 
to  include  a  smaller  point  target,  the  “knee”  of  the  curve  would  shift  to  a  smaller  CEP. 
Then  a  more  expensive,  0.5  m  CEP  weapon  might  become  more  cost  effective  then  less 
expensive,  1  m  CEP  rockets.  Conversely,  if  the  centered,  unbiased  aimpoint  assumption 
(discussed  on  page  10)  were  not  valid,  the  general  cost  effectiveness  advantage  of  the  high 
precision  rocket  would  become  smaller,  and  guided  rockets  with  slightly  larger  dispersion 
could  become  more  cost  efficient. 

Ideally,  COEA  studies  will  include  sensitivity  analyses  to  determine  the  range  of 
results  depending  on  changes  in  the  assumed  inputs.  With  these  sensitivities  in  mind,  the 
designer  can  use  the  COEA  models  as  powerful  design  tools. 
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in.  PRELIMINARY  SIZING  ANALYSES 


Several  aspects  of  any  design  require  preliminary  sizing.  For  the  NFS  HRGK 
design,  the  physical  dimensions  (length,  weight,  and  center  of  gravity),  the  thruster,  and 
the  seeker’s  field  of  view  each  needed  to  be  sized  early  in  the  design  process. 

The  author  developed  several  computer  codes  to  allow  the  sizing  of  the  kit,  its 
control  mechanisms,  and  the  seeker’s  field  of  view.  These  codes  each  use  top-level 
requirements  to  generate  lower-level  design  specifications.  This  flow  down  of 
specifications  helps  to  assure  that  the  final  design  meets  the  system  requirements. 

The  following  sections  describe  the  methodology  used  in  developing  the  computer 
models  and  describe  some  results  specific  to  the  NFS  HRGK  design.  The  chapter  ends 
with  observations  regarding  the  models  and  the  flow-down  of  design  requirements. 

A.  LENGTH,  WEIGHT,  AND  CENTER  OF  GRAVITY 

The  HRGK  design  weight  and  center  of  gravity  affects  the  overall  center  of  gravity 
of  the  guided  rocket.  The  guided  rocket’s  center  of  gravity  is  a  critical  design  consideration 
in  that  it  affects  the  rocket’s  stabihty  in  flight.  Additionally,  the  guided  rocket’s  center  of 
gravity  must  be  such  that  a  fully  loaded  launcher  pod’s  center  of  gravity  is  kept  between  the 
pod’s  two  mounting  lugs.  Furthermore,  the  length  and  weight  of  the  HRGK  are 
constrained  by  limits  specified  for  the  overall  length  and  weight  of  the  guided  rocket. 

A  Matlab  code  was  written  by  the  author  to  help  visualize  the  relationships 
between  the  HRGK  size  and  the  rocket’s  center  of  gravity.  The  code  also  maps  the  useable 
design  space  bounded  by  the  maximum  length  and  weight  as  well  as  the  launcher’s  center 
of  gravity  constraints.  This  code  is  briefly  described  and  some  results  are  presented  in  the 
following  subsections.  The  code  listing  is  included  in  Appendix  C. 

1.  Methodology 

The  computer  code  used  for  the  length,  weight,  and  center  of  gravity  analysis  is 
based  on  the  simple  equation  defining  a  composite  center  of  gravity,  X^q,  based  on  the 
component  weights,  W,,  and  component  centers  of  gravity,  x^g.  .  Namely, 

■  Iw,  ■ 
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a.  Inputs 

The  user-defined  parameters  in  the  code  are  either  design  or  constraint 
parameters  as  listed  in  Table  6.  Additional  model  parameters  needed  for  the  model  include 
the  weight  and  center  of  gravity  locations  for  the  live  and  fired  rocket  and  for  the  empty 
launcher.  The  length  of  the  rocket  and  the  location  of  the  pod’s  front  lug  are  also  required. 
These  parameters  for  the  Hydra-70  rocket  system  were  available  from  [Ref.  10  and  11]. 


Table  6.  Weight,  Length,  and  Center  of  Gravity  Model  Inputs. 


Input  Name 

Input  Description 

Design  Space  Parameter 

Wkl 

Range  of  kit  weights  (evenly-spaced  vector)  (lb) 

Lkl 

Specific  kit  lengths  to  be  plotted  (vector)  (in) 

Dkl 

Specific  kit  densities  to  be  plotted  (vector)  (Ibm/in^) 

Cgkl 

Assumed  kit  center  of  gravity  (fraction  of  kit  length  back  from  nose) 

dia 

Assumed  kit  diameter  (in) 

Constraint  Parameters 


Lmax 

Maximum  total  length  (in) 

Wmax 

Maximum  total  weight  (lb) 

CgL 

Pod  forward  Cg  limit  (fraction  of  lug  spacing  back  from  front  lug) 

The  kit’s  length  is  defined  as  the  length  extending  forward  from  the  tip  of 
the  rocket’s  fuze  to  the  forward  tip  of  the  kit.  The  kit’s  center  of  gravity  is  based  on  that 
length;  therefore,  if  the  kit’s  center  of  gravity  were  behind  the  tip  of  the  fuze,  the  value  for 
Cgkl  would  be  greater  than  unity. 

b.  Outputs 

The  computer  code  creates  a  plot  of  the  rocket’s  center  of  gravity’  as  a 
function  of  guidance  kit  weight.  The  design  space  is  marked  with  various  curves  indicating 
user  specified  guidance  kit  densities  and  lengths.  With  any  two  parameters  defined  (kit 
weight,  density,  length,  or  guided  rocket  center  of  gravity),  the  other  two  parameters  are 
fixed  on  the  plot.  The  usable  design  space  is  bounded  by  the  maximum  total  length  and 
weight  as  well  as  the  loaded  launcher’s  center  of  gravity  constraints.  The  code  generates  the 
described  plots  for  both  the  live  and  the  fired  rocket  cases. 


*  Measured  in  calibers  (missile  diameters)  from  the  forward  tip  of  the  guidance  kit. 
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2.  Results 

Figure  9  shows  the  output  from  the  computer  model  for  the  live  rocket  case.  The 
solid  lines  represent  the  case  where  the  kit’s  center  of  gravity  is  located  at  the  kit’s  defined 
mid-length.  The  cross-hatched  markings  around  the  lines  allow  for  the  kit’s  center  of 
gravity  to  travel  within  the  middle  third  of  the  kit’s  defined  length.  The  upper  edge  of  the 
cross-hatched  regions  represents  the  kit’s  center  of  gravity  being  two-thirds  of  the  length 
back  from  the  nose,  and  the  lower  edge  is  for  the  center  of  gravity  at  the  one-third  mark. 
The  bold  lines  represent  the  indicated  design  constraints.  The  NFS  HRGK  design  point  (7 
pounds  and  13.625  inches  longer  than  the  fuze)  is  marked  on  the  plot  as  a  reference. 


Figure  9.  Physical  Size  and  Center  of  Gravity  Design  Space.  Kit  length  is  defined  from  the 
tip  of  the  rocket  fiize  forward.  The  launcher  center  of  gravity  constraint  is  tighter  than 
specified  for  illustrative  purposes. 


From  the  plot  it  appears  that  the  loaded  launcher’s  center  of  gravity  requirement 
(keeping  the  pod’s  center  of  gravity  behind  the  launcher’s  front  mounting  lug)  is  not 
binding  for  most  conditions.  Indeed,  for  Figure  9,  the  constraint  was  tightened  to  keeping 
the  pod’s  center  of  gravity  within  the  middle  third  of  the  14  inch  lug  spacing  so  that  the 
constraint  could  even  be  seen  in  the  plot. 
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B.  MANEUVERABILITY 


Another  computer  model  was  developed  in  Matlab  to  define — based  on  the  system 
requirements — ^the  HRGK’s  maximum  maneuverability  specifications.  An  additional 
computer  program  uses  a  steady  maneuver  to  determine  the  angles  of  attack  and  thruster 
forces  needed  to  achieve  the  specified  maneuverability.  This  latter  code  sizes  the  HRGK’s 
thrusters.  The  following  subsections  describe  these  computer  models  and  their  results. 

1.  Turn  Rate  and  Turn  Radius  Methodology 

The  maneuverability  analysis  determines  the  maximum  turn  rate  in  g’s  or  minimum 
turn  radius  in  kilometers  needed  to  hit  a  target  over  various  launch  and  seeker  acquisition 
ranges.  The  two-dimensional  model  used  in  this  analysis  is  described  in  this  subsection. 

a.  Model  Inputs 

The  maneuverability  model  requires  several  user-defined  inputs.  These 
inputs  are  listed  and  described  in  Table  7  and  are  illustrated  graphically  in  Figure  10. 


Table  7.  Maneuverability  Model  Inputs. 


Input  Name 

Input  Description 

Launch  and  Seeker  Acquisition  Ranges 

Rmin,  Rmax 

Minimum  and  maximum  range  to  target  at  launch  (Ion) 

r% 

•*^acq 

Specific  seeker’s  maximum  target  acquisition  range  (vector)  (km) 

Unguided  Errors 

Pnt 

Launcher  pointing  error  at  time  of  launch — la  (degrees) 

B 

Ballistic  unguided  flyout  error — la  (milliradians) 

Weapon  System  Delays  and  Limitations 

Gd 

Guidance  delay  (time  after  launch  before  guidance  is  possible)  (sec) 

Target  identification  delay  (from  initial  acquisition  until  start  of  turn)  (sec) 

Sp 

Straight  final  flight  distance  (terminal  non-maneuvering  zone)  (m) 

Missile  and  Target  Motion 

Time 

Time  values  for  missile  range  and  velocity  time  profiles  (vector)  (sec) 

Dist 

Missile  range  profile  (distances  corresponding  to  Time)  (vector)  (m) 

Vm 

Missile  velocity  profile  (velocities  corresponding  to  Time)  (vector)  (m/s) 

Vt 

Target  speed  (mph) 

Statistical  Confidence 

sig 

1 

Number  of  sigma’s  confidence  (e.g.,  sig=2  gives  ±2a  or  95%  confidence) 
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Figure  10.  Graphical  Representation  of  Maneuverability  Model  Inputs. 


The  two-dimensional  model  assumes  that  the  rocket  flies  a  straight  path 
toward  the  initial  target  location  but  with  an  azimuth  error  angle,  6.  This  angle  is  a 
statistical  combination’  of  the  launcher  pointing  error  and  the  ballistic  flyout  dispersion 
error  of  the  rocket.  After  both  the  initial  guidance  delay^,  Gp,  and  after  the  rocket  is  within 
the  assumed  seeker  acquisition  range,  of  the  target;  the  rocket  can  begin  a  turn  toward 

the  target — correcting  any  cross-track  errors  due  to  6ot  target  motion. 

A  target  identification  delay  time,  tp,  between  the  first  detection  of  the  target 
and  the  start  of  the  turn  can  also  be  included  in  the  model.  The  distance  associated  with  tp, 
Sp,  is  based  on  the  rocket  velocity  at  the  time  of  first  detection  (or  the  end  of  guidance 
delay,  whichever  is  later).  The  model  also  allows  for  a  terminal  non-maneuver  zone — some 
radius,  Sp  from  the  target.  For  modeling  purposes,  this  zone  around  the  target  forces  the 
rocket  to  have  completed  its  turn  to  the  target  and  align  itself  for  a  final  straight  approach. 

b.  Stationary  Target — Closed-Form  Solution 

With  a  stationary  target,  the  inputs  described  above  can  provide  a  closed- 
form  solution  for  the  maneuverability  specifications.  This  process  is  laid  out  in  Figure  11. 

The  equations  illustrated  in  Figure  1 1  provide  the  maneuverability 
requirements  in  terms  of  both  required  turn  radius,  r^,  and  turn  rate,  n^.  (They  also 
provide  the  solution  to  the  minimum  required  seeker  field  of  view — FOV.  This  part  of  the 
analysis  will  be  described  more  fully  in  a  later  section.) 

For  the  ±2cr  unguided  error  case  {sig  =  2),  the  resulting  maneuverability 
values  (r^-and  n^)  could  be  interpreted  as  the  biggest  turn  radius  or  slowest  turn  rate  that 
would  still  ensure  a  95%  probability  of  hitting  a  stationary  target.  In  other  words,  any 
larger  or  slower  turns,  than  those  specified,  would  not  give  a  95%  confidence  with  the 
provided  inputs.  On  the  other  hand,  tighter  or  faster  turns  would  give  a  higher  than  95% 
confidence.  Therefore,  if  the  system  is  specified  to  95%  confidence,  the  computed  values 
are  the  maximum  turn  rate  and  minimum  turn  radius  required. 

c.  Moving  Target 

The  basic  equations  shown  in  Figure  1 1  can  be  applied  to  the  moving  target 


'  The  standard  deviation,  a,  of  6  is  the  square  root  of  the  sum  of  the  variances  Pnf  and  — 
I  2  2 

a  =  iPnt  +B  — ^where  Pnt  and  B  are  the  standard  deviations  for  the  pointing  error  and  the  ballistic 
dispersion  errors,  respectively.  The  angle  6  is  some  user-defined  multiple  (sig)  of  CT. 

^  The  range  associated  with  is  calculated  based  on  the  user-supplied  range  time  profile  for  the  rocket. 
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Figure  11.  Closed-Form  Solution  for  Maneuverability  Specifications  with  a  Stationary 

Target. 
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scenario  shown  in  Figure  12.  The  moving  target  scenario  assumes  a  “worst-case”  target 
motion.  The  target  moves  with  a  crossing  velocity  in  the  direction  opposite  the  flyout  error. 
Then,  after  it  is  detected,  the  target  turns  in  at  45°  making  the  rocket’s  turn  tighter. 

The  moving  target  scenario  precludes  a  closed-form  solution  of  the 
maneuverability  specifications.  However,  the  equations  in  Figure  1 1  can  be  used  with  an 
iterative  step  at  each  of  the  flight  milestones  (target  detection,  start  of  turn,  and  end  of  turn). 
The  Matlab  code  listing  in  Appendix  C  uses  this  approach.  The  following  description 
provides  an  example  of  the  iterative  process  used  in  the  computer  code. 

After  Rp  is  calculated  using  the  appropriate  stationary  target  equations  (see 
Figure  11),  the  time  required  to  fly  that  distance  is  available  from  the  range  time  profile. 
Using  the  time  and  the  target  speed,  the  new  target  location  can  be  computed.  The  value  for 
Rd  can  then  be  re-calculated  based  on  the  updated  target  location.  These  iterations  continue 
until  the  change  in  between  iterations  is  less  than  some  set  tolerance  (such  as  1  m). 

2.  NFS  HRGK  Turn  Rate  and  Turn  Radius  Analysis 

This  subsection  discusses  the  maneuverability  analysis  performed  specifically  for 
the  design  of  the  NFS  HRGK.  The  following  paragraphs  describe  the  inputs  used  for  the 
computer  model  and  present  the  analysis  results. 

a.  NFS  HRGK  Model  Inputs 

The  NFS  HRGK  analysis  used  the  full  gamut  of  specified  launch  ranges  (1 
to  6  kilometers)  and  various  seeker  acquisition  ranges  (0.75,  1.0,  1.15,  2.0,  3.0,  and  6.0 
km).  The  assumed  la  values  for  launcher  pointing  error  and  ballistic  dispersion  were  y 
and  12  mils,  respectively.  (Using  a  95%  confidence  level — sig  =  2,  the  angular  flyout  error 
was  internally  computed  to  be  6.16°.)  The  target  speed  was  60  mph  per  the  system 
requirements.  Initial  guidance  was  delayed  as  a  linear  function  of  launch  range  (0.5 
seconds  for  1  km,  up  to  1.0  second  for  6  km).  The  target  identification  delay  was  0.625 
seconds  (arbitrarily  chosen  to  correspond  to  5  laser  pulses  at  8  Hertz),  and  the  terminal 
non-maneuvering  distance  was  arbitrarily  chosen  to  be  2(X)  meters. 

The  rocket’s  velocity  and  distance  time  profiles  used  in  the  analysis  were 
based  on  an  early  simulation  of  a  HRGK-equipped  rocket  using  a  ballistic  trajectory. 
Specifically,  the  simulation  used  the  150°F  propellant  thrust  profile  [Ref.  1 1]  and  a  generic, 
hemispherical-nosed  guidance  kit  with  three  canards.  The  aerodynamic  coefficients  and 
derivatives  used  in  the  simulation  were  determined  using  the  software  code.  Advanced 
Design  of  Aerodynamic  Missiles — ^ADAM  [Ref.  12] . 
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Figure  12.  Moving  Target  Scenario  for  Maneuverability  Analysis. 


b.  NPS  HRGK  Analysis  Results 

Based  on  the  inputs  described  above,  the  maximum  turn  rate  and  minimum 
turn  radius  specifications  for  the  NPS  HRGK  are  shown  in  Figure  13  and  Figure  14, 
respectively.  The  results  are  shown  as  a  function  of  launch  range  for  several  seeker 
acquisition  ranges.  The  analysis  results  show  that  as  long  as  the  seeker  can  initially  detect 
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the  target  at  1.15  kilometers  or  more,  the  hardest  g-loading  maneuver  for  the  rocket  would 
be  less  than  15  g’s  and  tightest  turn  radius  would  be  greater  than  1.83  kilometers. 

The  maximum  turn  rate  is  a  strong  function  of  launch  range  with  short  range 
launches  requiring  higher  g-loading  maneuvers.  This  is  primarily  due  to  higher  rocket 
speeds  at  shorter  ranges.  For  example,  after  flying  0.5  km,  the  rocket’s  velocity  is  1.8 
times  higher  than  it  is  after  coasting  out  to  5  km.  For  equal  turn  radii,  a  1.8  factor  increase 
in  velocity  yields  a  3.2  (or  1.8^)  factor  increase  in  turn  rate.  The  guidance  delay  also 
increases  the  turn  rate  at  shorter  launch  ranges.  Target  identification  delays  primarily  impact 
the  shorter  seeker  acquisition  range  cases — longer  delays  causing  higher  turn  rates. 

Unlike  the  maximum  required  turn  rate,  the  minimum  required  turn  radius  is 
almost  independent  of  rocket  speed'.  Once  the  launch  range  places  the  target  initially 
outside  the  seeker’s  acquisition  range,  the  turn  radius  curves  remain  fairly  flat  with 
increasing  launch  range.  For  example,  with  a  1.15-km-acquisition-range  seeker,  the 
minimum  required  turn  radius  changes  less  than  59%  across  the  spectrum  of  launch  ranges 
(1.83  km  at  1  km  launch  to  2.90  km  at  6  km  launch).  Conversely,  if  the  seeker  has 
sufficient  detection  range  to  take  advantage  of  it,  the  longer  range  launches  allow 
significantly  wider  turns.  For  example,  only  a  gentle  12.9  km  turn  radius  is  required  with  a 
6  km  seeker  acquisition  range  and  a  6  km  launch  range. 

By  either  measure  (turn  rate  or  turn  radius),  the  most  stressing  case  for  the 
maneuverability  requirement  is  the  short  range  launch,  because  of  rocket  speed  and  the  lack 
of  reaction  time.  The  seeker’s  acquisition  range  can  also  significantly  impact  the  HRGK 
maneuverability  requirements — shorter  range  seekers  forcing  higher  maneuverability 
requirements.  Conversely,  long  range  seekers  can  alleviate  some  maneuverability 
requirements  by  detecting  the  target  earlier  and  allowing  gentler  turns. 

An  additional  benefit  of  the  maneuverability  analysis  is  that  the  minimum 
required  seeker  acquisition  range  is  shown  to  be  1.15  km.  With  shorter  acquisition  ranges 
and  long  launch  ranges,  the  target’s  cross-track  position  relative  to  the  rocket  (due  to  0  and 
target  motion)  exceeds  the  seeker’s  range.  For  example,  with  a  4.1  km  launch  range,  the 
rocket’s  closest  unguided  approach  to  the  target  is  0.75  km  (2a).  For  longer  launch  ranges, 
the  closest  approach  distance  would  be  greater  than  0.75  km;  therefore,  a  0.75-km- 
detection-range  seeker  would  not  have  a  95%  probability  of  detecting  the  target. 

'  Rocket  speed  indirectly  affects  the  turn  radius  since  it  partially  determines  the  distance  the  target  can 
travel  before  it  is  detected  or  impacted. 
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Figure  13.  Maximum  Turn  Rate  Required  to  Hit  Target  (95%  Probability).  Only  with  a 
seeker  able  to  detect  a  target  from  1.15  km  or  more  >  1.15)  can  the  maximum 
range  requirement  (6  km)  be  meet. 
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_ _ _ Launch  Range  (km) 

Figure  14.  Minimum  Turn  Radius  Required  to  Hit  Target  (95%  Probability).  Seeker 
acquisition  ranges  are  indicated  with  the  corresponding  curves. 
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3.  Steady-Turn  Maneuver  Analysis 

Once  maneuverability  specifications  are  set,  a  steady-tum  maneuver  model  can  be 
used  to  size  the  control  mechanism  needed  to  achieve  the  needed  turn  rate  or  turn  radius. 
The  following  paragraphs  describe  a  computer  code  developed  by  the  author  to  perform 
this  task.  Later  a  specific  example  for  a  guided  rocket  is  presented.  The  Matlab  code 
listing  with  the  NFS  HRGK  parameters  is  included  in  Appendix  C. 


a.  Mathematical  Steady-Turn  Model 

The  purpose  of  this  mathematical  steady-tum  model  is  to  determine  the 
thrast  force  and  angle  of  attack  (or  sideslip  angle)  necessary  to  generate  a  constant 
acceleration  turn  in  the  horizontal.  Throughout  the  description,  the  more  familiar 
longitudinal  nomenclature  (loosely  associated  with  the  vertical  plane)  is  used  even  though 
the  turn  is  described  in  the  directional  case  (loosely  associated  with  the  horizontal  plane). 
The  interchange  of  nomenclature  is  justified  since,  for  the  guided  rocket  or  missile  case,  the 
longitudinal  and  directional  dynamics  are  equivalent.  A  useful  listing  of  missile 
aerodynamic  nomenclature  is  included  in  Appendix  A. 

Figure  15  shows  a  free-body  diagram  for  a  horizontal  maneuver.  For  a 
constant  acceleration  turn,  the  sum  of  the  forces  and  moments  can  be  expressed  by  the 
following  equations: 


Forces  =  nW  =  N  -i-  T  =  +  kT 


2^0 


pV^AC. 


a+kT 


J 


^Moments  =  -l- TXj.  - =  0 

=  —{^qACf]^c^X(^p  +  kTXj  —  qAc^C/^^  ^ 


-  -  - 


’  2V 


=  -[qACff  C^X(,p  +  kTXp  -  qAd 
=  -{qAC„a)x^  +  *7X, 


=  a\ 


\Xcp  +  TkXp  --^pACu  d^ng 
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where  ^  is  a  thrust  multiplication  factor  (accounting  for  misalignment  of  the  thruster  and  the 
effect  of  jet  and  free  air  steam  interaction)  and  the  other  terms  are  either  defined  by  Figure 
15  or  are  standard  terms  defined  in  Appendix  A.  The  moments  equation  uses  the  fact  that 


rj. 


V 

) 


Figure  15.  Plane-View,  Free-Body  Diagram  for  a  Steady-Turn  Maneuver  (horizontal,  left 
turn — gravity  force  vector  pointed  into  the  page). 


The  sum  of  forces  and  sum  of  moments  equations  can  be  rewritten  as  the 
following  matrix  equation: 


^8o 

'T' 

nW 

(  1 

-—PV^AC,^  X,, 
1  )  \ 

a 

-^pAC^d^ng 

l^So  ’  J 

which  gives  the  thrust  magnitude  and  angle  of  attack  in  terms  of  guided  rocket  parameters 
and  the  turn  rate  n.  Using  the  relationship  =  VVitg,  the  matrix  equation  can  be  expressed 
in  terms  of  the  turn  radius  as  the  following: 
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b.  Steady-Turn  Computer  Code 

The  author  has  developed  a  computer  code  to  implement  the  mathematical 
model  described  above.  This  code  solves  for  the  thrust  and  angle  of  attack  required  to 
achieve  the  required  turn  rate  or  turn  radius  over  a  range  of  Mach  numbers.  The  code 
requires  the  user-supplied  inputs  listed  and  described  in  Table  8. 


Table  8.  Steady-Turn  Maneuver  Computer  Model  Inputs. 


Input  Name 

Input  Description 

Scenario  Definition 

n  or 

Desired  turn  rate  (g’s)  or  turn  radius  (ft) 

type 

Indicator  for  computations  based  either  on  turn  rate  (1)  or  turn  radius  (2) 

alt 

Pressure  altitude  (ft) 

Mmin,  Mmax 

Minimum  and  maximum  Mach  number  for  the  calculations 

Rocket  Configuration 

Wght 

Rocket  total  weight  at  maneuvering  time  (lb) 

Dia 

Rocket  diameter  (for  reference  area  and  length)  (in) 

Xcg 

Rocket  center  of  gravity  (calibers  from  nose) 

Xt 

Thruster  location  (calibers  from  nose) 

Aerodynamic  Parameters  (functions  of  Mach) 


Mach_dat 

Mach  values  for  corresponding  aerodynamic  parameters  (vector) 

CNa_dat 

Rocket  normal  force  curve  slope  (corresponding  to  Mach)  (vector)  (rad  *) 

Xcp_dat 

Center  of  pressure  (corresponding  to  Mach)  (vector)  (calibers  from  nose) 

CMq_dat 

“Pitch”  damping  (corresponding  to  Mach)  (vector)  (rad"') 

c. 

NFS  HRGK  Example 

The  following  provides  example  solutions  to  the  equations  previously 
developed  for  the  steady-tum  model.  The  examples  are  specific  to  the  NFS  HRGK  with  the 
guided  rocket’s  maneuverability  specifications  set  at  15  g’s  turn  rate  or  1.83  km  (6,000  ft) 
turn  radius.  The  manual  calculations  are  for  the  NFS  HRGK  design  at  Mach  1.0.  These 
calculations  are  done  first  for  the  turn  rate  specification  and  then  for  the  turn  radius 
specification.  The  aerodynamic  coefficients  are  from  ADAM  [Ref.  12].  All  computations 
are  for  500  feet  pressure  altitude  and  use  the  fired  rocket  center  of  gravity  (most 
maneuvering  takes  place  after  the  1.0  sec  rocket  motor  bum).  Because  no  additional  data  is 
available,  the  value  of  k  (an  adjustment  factor  for  the  thruster  force)  is  assumed  to  be  unity. 
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The  manual  calculations  are  followed  by  sample  output  from  the  steady-tum  computer 
code. 

( 1 )  Sample  Manual  Calculation  Based  on  Maximum  Turn  Rate. 


1 

(0.0754f^)(l  1 14f )'(0.0412  ft')(20.05  rad  ' ) 

'T 

2(32.2  "’'“•M 

V  Ibf  s^  / 

1.79ft 

(0.0754^)(1 1 14f )'(0.0412  ft")(20.05  rad  ') ^ 

/  „  \  0.0824  ft 

0/  ao  0  Ibrn-ft  \ 

Ibf-s^  J 

a 

15  (22.8  Ibf) 

(0.0754^)(0.0412  ft')(-5520  rad-')(0.229  ft)'l5  (32.2^) 
The  solution  to  the  set  of  equations  is: 


1 

1200M- 

'T' 

■  342  Ibf  ■ 

1.79  ft 

-99.1^_ 

a 

-3.38  Ibf 

T 

'  13.2  Ibf  ■ 

‘13.2  Ibf 

a 

0.273  rad 

.  15.6°  _ 

(2)  Sample  Manual  Calculation  Based  on  Minimum  Turn  Radius. 


(I114ff  (22.8  Ibf) 

(32.2^)(6000  ft) 

(0.0754-|^)(lll4f)^ (0.0412  ft^)(-5520  rad“’)(0.229  ft)^ 
4(32.2  ^)(6000  ft) 

The  solution  to  the  set  of  equations  is: 


1200^ 


1.79  ft  -99.11^ 

rad  . 


a 


1  1200^  ■ 

'T' 

■  147  Ibf  ■ 

1.79  ft  -99.1^ 

a 

-1.45  Ibf -ft 

T 

■  5.68  Ibf  ■ 

■5.68  Ibf 

a 

0.117  rad 

6.71° 

37 


(3)  Sample  Output  from  Steady-Tum  Computer  Code.  Figure 
16  shows  the  output  from  a  computer  code  developed  by  the  author.  The  code  computes 
the  solution  to  the  steady-tum  equations  that  were  manually  solved  in  the  previous 
paragraphs. 


Figure  16.  Required  Thrust  and  Angle  of  Attack  for  Maximum  Turn  Rate  (15  g's — 
dashed  line)  or  Minimum  Turn  Radius  (1.83  km — solid  line)  (NFS  HRGK  design). 


From  the  data  plotted  in  Figure  16,  the  difference  in  the  two 
methods  of  specifying  maneuverability  requirements  is  apparent.  For  example,  at  subsonic 
speeds,  the  model  predicts  that  the  rocket  could  only  achieve  a  15  g  maneuver  with  very 
high  thrust  and  extremely  high  angle  of  attack.  On  the  other  hand,  the  requisite  1.83  km 
turn  radius  can  be  achieved  at  subsonic  speeds  with  a  moderate  angle  of  attack  and  low 
thruster  force. 

Fortunately,  15  g  maneuvers  are  not  required  when  the  missile  has 
coasted  down  to  subsonic  speeds  (see  long  range  side  of  Figure  13,  page  33).  Clearly,  the 
reasonable  curve  for  sizing  the  thruster  is  the  minimum  turn  radius  curve  (solid  line).  This 
seems  reasonable  since  the  turn  radius  requirements  remained  fairly  flat  over  the  entire 
gamut  of  launch  ranges  (see  Figure  14). 
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Designing  to  the  minimum  turn  radius  and  given  the  fact  that  the 
rocket  has  slowed  down  to  nearly  Mach  1.3  by  the  first  kilometer  of  flight,  the  thruster  can 
be  sized  to  provide  10  to  15  pounds  of  thrust.  A  design  with  higher  thrust  early  in  the  flight 
would  be  ideal,  since  the  short-range,  high-speed  scenario  is  the  most  demanding. 

4.  Probabilistic  Approach 

The  maneuverability  analysis  performed  in  the  previous  subsections  was  based  on 
the  statistical  (probabilistic)  flyout  of  the  rocket  but  a  deterministic,  worst-case  target 
motion.  Another  approach  to  determining  the  needed  level  of  maneuverability  would  be  to 
apply  a  statistical  model  to  the  target  motion  as  well  as  the  flyout.  This  latter  approach 
would  give  valuable  insight  into  the  reasonableness  of  the  worst-case  approach  taken 
above.  This  type  of  purely  probabilistic  approach  was  needed  for  the  seeker  field  of  view 
analysis,  which  is  described  next. 

C.  SEEKER  FIELD  OF  VIEW 

The  seeker’s  field  of  view'  was  another  early  consideration  for  the  NPS  HRGK 
design.  Again,  this  design  specification  could  be  derived  fi-om  top-level  system 
requirements  and  simple  mathematical  models.  A  possible  approach  to  determining  the 
needed  seeker  field  of  view  was  alluded  to  in  the  maneuverability  analysis  section.  This 
section  will  briefly  discuss  that  approach  and  its  results.  Then  a  probabilistic  approach  with 
its  results  will  be  discussed  in  greater  detail. 

1.  Analysis  with  “Worst-Case”  Target  Motion  Model 

The  two-dimensional  scenario  described  in  subsection  in.B.l.  Turn  Rate  and  Turn 
Radius  Methodology  provides  a  means  for  calculating  the  needed  HRGK  seeker  field  of 
view.  The  equations  shown  in  Figure  11  (page  29),  can  be  used  to  determine  the  field  of 
view  that  keeps  a  stationary  target  within  the  seeker’s  field  of  view  with  some  specified 
level  of  confidence  (for  example,  95%)  up  until  the  target  can  be  detected  and  a  turn  toward 
the  target  can  be  initiated.  As  discussed  in  the  earher  subsection,  the  equations  shown  in 
Figure  1 1  can  be  solved  iteratively  to  determine  the  required  field  of  view  for  a  moving 
target  (see  Figure  12,  page  31). 


‘  The  seeker’s  field  of  regard— meaning  the  total  area  covered  by  a  fixed  or  scanning  seeker— is  the  actual 
consideration,  but  for  a  strapped-down  seeker,  the  field  of  view  and  field  of  regard  can  be  used 
interchangeably.  This  thesis  uses  the  term  field  of  view. 
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Unfortunately,  the  target  motion  model  used  for  the  maneuverability  analysis  is  a 
worst-case  example.  (The  target  crosses  at  maximum  velocity  in  the  opposite  direction  of 
the  statistical  flyout  and  pointing  errors.  Furthermore,  after  detection,  the  target  turns  in  at 
45°  which  creates  the  need  for  a  still  larger  field  of  view.  See  Figure  12.) 

The  required  field  of  view  results  shown  in  Figure  17  for  the  NFS  HRGK  were 
obtained  using  the  worst-case  target  motion  model  and  the  computer  code  inputs  discussed 
in  subsection  III.B.2.  NPS  HRGK  Turn  Rate  and  Turn  Radius  Analysis  (page  30).  The 
values  are  for  half  field  of  view  and  are  shown  as  a  function  of  launch  range  and  for 
various  seeker  acquisition  ranges.  For  example,  with  a  maximum  launch  range  and  a  seeker 
acquisition  range  of  1.5  km,  a  ±48°  seeker  field  of  view  would  be  needed  to  assure  (95%) 
that  the  target  stays  in  the  field  of  view  up  until  a  turn  can  be  initiated.  (With  a  stationary 
target,  only  a  ±29°  seeker  field  of  view  would  be  needed.) 


Figure  17.  Required  Field  of  View  Results  for  95%  Probability  of  Target  Staying  in  Field 
of  View  with  Worst  Case  Target  Motion.  (Corresponding  field  of  view  requirements  for  a 
stationary  target  are  indicated  with  dashed  lines.) 

From  the  curves  in  Figure  17,  it  is  clear  that  longer  launch  ranges  require  larger 
fields  of  view  (after  longer  flyouts,  the  rocket’s  cross-track  error  and  the  target’s  motion 
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away  from  the  rocket  will  be  larger).  Shorter  acquisition  range  seekers  also  require  larger 
fields  of  view. 

Unfortunately,  the  indicated  field  of  view  sizes  for  long  range  launches  and  short 
range  seekers  are  unattainably  large,  especially  for  a  strapped-down,  semi-active  laser  spot 
homing  seeker  like  that  on  the  NPS  HRGK.  The  worst-case  target  motion  model  is 
unacceptable,  and  so,  a  probabilistic  approach  must  be  taken. 

2.  Analysis  with  Probabilistic  Target  Motion  Model 

The  probabilistic  target  motion  model  used  for  the  seeker  field  of  view  analysis  is 
based  on  a  distribution  of  target  positions.  In  this  analysis,  the  probability  of  the  seeker 
acquiring  a  target  is  derived  from  an  assumed  seeker  field  of  view  and  the  probabilistic 
distributions  for  the  rocket  and  the  target  The  methodology  for  defining  a  seeker  field  of 
view  specification  with  this  approach  is  developed  in  this  subsection.  The  resulting 
specifications  for  the  NPS  HRGK  are  also  presented.  The  Matlab  code  developed  by  the 
author  to  perform  this  analysis  is  included  in  Appendix  C. 

a.  Computer  Model  Inputs 

The  inputs  required  for  the  probabilistic  seeker  field  of  view  analysis  are 
listed  and  described  in  Table  9. 


Table  9.  Seeker  Field  of  View  Analysis  Inputs. 


Input  Name 

Input  Description 

Seeker  Descriptions 

FOVmin,  max 

Minimum  and  maximum  seeker  half  field  of  views  (deg) 

^acq 

Specific  seeker’s  maximum  target  acquisition  range  (vector)  (km) 

Unguided  Errors  and  Weapon  System  Delays 


Pnt 

Launcher  pointing  error  at  time  of  launch — lo  (degrees) 

B 

Ballistic  unguided  flyout  error — lo  (milliradians) 

Go 

Guidance  delay  (time  after  launch  before  guidance  is  possible)  (sec) 

Target  identification  delay  (from  initial  acquisition  until  start  of  turn)  (sec) 

Missile  and  Target  Motion 


RL 

Range  to  target  at  launch  (single  value)  (km) 

Time 

Time  values  for  missile  range  and  velocity  time  profiles  (vector)  (sec) 

Missile  range  profile  (distances  corresponding  to  Time)  (vector)  (m) 

^Tmin»  ^Tmax 

Minimum  and  maximum  target  speed  (mph) 
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b.  Target  Motion  Model 

The  probabilistic  seeker  field  of  view  analysis  uses  a  simple  distribution 
model  for  target  motion.  The  target’s  speed  is  assumed  to  be  constant  and  uniformly 
distributed  between  the  user-supplied  inputs,  and  V,.^.  Its  bearing,  (j),  is  also 

assumed  to  be  constant  and  to  be  uniformly  distributed  between  (f  and  360°.  Figure  18 
shows  a  plane  view  of  possible  target  locations  after  time,  tj  (the  flyout  time  before  the 
rocket  starts  its  turn  to  the  target).  The  figure  also  includes  the  corresponding  probability 
density  function  for  the  target’s  cross-track  position,  jc.  This  function  provides  the 
probability  of  the  target  being  at  some  specific  cross-track  position,  Xj,  or  between  two 
cross-track  limits,  Xj  and  X2.  (The  y-axis  indicates  the  direction  of  the  rocket  flyout.) 


Figure  18.  Target  Distribution  Model,  (a)  Plane  view  of  possible  target  positions  and  (b) 
Cross-track  probability  density  function,  Pi(x) 


c.  Analysis  Methodology 

The  methodology  used  to  compute  the  probability  of  the  target  being  within 
the  HRGK  field  of  view,  has  previously  been  documented  by  the  author.  The  following 
paragraphs  briefly  summarize  the  methodology.  In  these  paragraphs,  the  assumption  is 
made  that  if  the  target  is  within  the  HRGK’s  field  of  view,  it  will  be  acquired;  therefore,  the 
desired  quantity  is  the  probability  of  acquisition,  The  referenced  technical 

memorandum  should  be  consulted  for  a  more  complete  description  of  the  theory  and  the 
limitations  of  the  simplifying  assumptions.  [Ref.  13] 

Figure  19  shows  a  sample  case  for  the  probabilistic  model  and  shows 
several  of  the  variables  used  in  the  model  computations.  To  obtain  the  probability  of  the 
target  being  within  the  seeker  field  of  view  {x,  <  <  x^,  the  independent  probabilities  of 
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the  rocket  and  target  cross-tracks  must  be  computed.  The  rocket  cross-track  computation 
assumes  a  normal  distribution  of  0 — which  depends  only  on  the  independent,  normally 
distributed  pointing  and  ballistic  flyout  errors,  Pnt  and  B,  and  the  target's  position  depends 
on  the  target  motion  model  described  above.  The  cross-track  distribution  for  the  rocket  is 
evaluated  for  the  assumed  range  to  detection  (/?o),  and  the  target  distribution  is  evaluated 
for  the  assumed  time  at  start  of  turn  {tj,  time  required  for  the  rocket  to  fly  to  Rj). 


Figure  19.  Illustrative  Case  for  Probabilistic  Seeker  Field  of  View  Analysis.  The  rocket  is 
shown  at  the  assumed  position  /?„. 


The  probability  of  the  rocket  being  at  cross-track  position  Xj^  is 


Pr{x)  = 


-x^ 

,2(72 


<J-J^ 


where  c  —  tan0  —  /?D-\/tan^(Pnf)-i-  .  Because  of  the  assumed  target  distribution,  the 

target  has  equal  probability  of  being  anywhere  in  the  possible  region  shown  in 

Figure  20.  Therefore,  the  probability  of  the  target  being  within  the  field  of  view  limits  {Xj 
and  X2)’ — ^the  shaded  area.  A,  in  the  figure — is  the  ratio  of  the  seeker  acquired  area,  A,  to 
the  total  possible  area,  A^qj.  The  seeker  acquired  area.  A,  can  be  expressed  as  the  area 
between  the  limits  and  x^  and  inside  the  outer  circle  (radius  =  minus  the  area 
between  the  limits  and  inside  the  inner  circle  (radius  =  Vj).  The  total  possible  area,  Ajqj,  is 
^{^2  “  ^i^)»  ti^us,  the  desired  probability  is  given  as: 


‘  The  field  of  view  limits  are  a:;  =  -  F  and  ;t2  =  +  F  where  the  offset,  F  =  F„,  sin(FOV). 
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j(gn(x2)min(^J*2])__  j/gn(x2)niin(/iJ*2|)__ 

2  j^|^2  dx  -  2  —x^dx 

p ^  \  _  «*n(*i)min(r2,|.>t,|)  iigrt(i,)min(r,,|j:,|) 

^OT  j 

where  the  limits  of  the  two  integrations  are  set  to  never  exceed  the  radius  of  the  outer  or 
inner  circle,  respectively. 


Figure  20.  Possible  Target  Locations.  The  ratio  of  the  area  (A)  bounded  by 
the  field  of  view  limits — x,  and  jcj — and  the  total  possible  area  (Aj-ot)  is  the 
probability  of  the  target  being  within  the  two  limits. 


An  alternative  calculation  of  the  seeker  acquired  area.  A,  is  illustrated  in 
Figure  21.  The  calculation  method  relies  on  the  fact  that  since  the  probabilities  are 
symetrical  about  the  y-axis,  the  calculations  only  need  to  be  made  for  positive  values  of  Xj,  ; 
thus  Xj  will  always  be  positive.  The  method  involves  calculating  the  area,  A,,  bounded  by 
the  circles,  the  lower  field  of  view  limit  (Xj),  and  the  y-axis  as  given  by  the  .equation: 

/i=min(r2,JC2) _  L2=imn(r|,jC2) _ 

Aj=2  jyjr^  -x^dx-2  -x^dx 

0  0 


=  AV''/-4+''2 


sm 


-1 


\ 


-  sm~^ 


fr  \ 


V'i  J 


Next,  the  area  Aj  is  calculated.  This  area  is  bounded  by  the  circles,  the  lower  field  of  view 
limit  (x,),  and  the  y-axis,  and  is  given  by  the  equation 
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/3  =min(r2  >|a:J)  L4  =nun(ri 

A2  =  2  -x^dx-  2  J-y/r,^  -jc^dlx 

0  0 


fr  \ 

=  I^^r^-L\  +r2^sin“‘ 

ZL 

-Ll  sin  ^ 

\ 

I'l  J 

Finally,  the  area  A2  is  subtracted  from  A,  if  aTj  is  greater  than  zero  (Figure  21  (a))  or  added 
to  Aj  ifx,  is  less  than  zero  (Figure  21  (b))— A  =A,  -  sign(A:,)  A2. 


Figure  21.  Alternative  Method  for  Computing  Target  Distribution  Area  within  Seeker 
Field  of  View  Limits.  (A  =  Ai  -  sign(x,)  Aj.) 


Now,  the  probability  of  the  rocket  being  at  a  specific  cross-track  and  the 
target  being  within  the  HRGK’s  field  of  view  can  be  expressed  as  the  joint  probability 
Xj<Xj<  X2)  =  Pji(Xji)  Pj(Xj,X2).  The  total  probability  of  the  rocket  having 
the  target  in  its  field  of  view  for  any  flyout  is  the  integral  of  this  joint  probability.  These 
probabilities  are  illustrated  in  Figure  22.  The  actual  method  for  computing  the  probability  of 
acquisition  is  described  in  the  following  paragraphs. 


Figure  22.  Probability  Density  Functions,  (a)  Rocket  cross-track  error,  (b)  Target  cross¬ 
track  position,  and  (c)  Probability  of  acquisition. 
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d.  Computer  Code  Description 

The  Matlab  computer  code  developed  to  determine  the  probability  of  target 
acquisition’  includes  a  Matlab  script  file  and  a  Matlab  function.  (Appendix  B  briefly 
describes  the  difference  between  script  files  and  functions.)  The  script  file  sets  the  user 
inputs,  determines  the  rocket’s  cross-track  error  standard  deviation,  sets  the  seeker  offset 
{F),  calculates  the  integral  of  the  joint  probability  \Pg{x)  P^x±F)dx  from  zero  to  an  upper 
bound,  and  plots  the  results.  The  Matlab  function  is  called  by  the  integration  scheme  in  the 
script  file  and  returns  values  for  the  joint  probability  =  P^{x)  P^x±F)).  These 

functions  are  all  performed  over  a  range  of  fields  of  view  and  for  several  seeker  acquisition 
ranges.  Both  the  script  file  and  the  function  are  included  in  Appendix  C. 

e,  NPS  HRGK  Results 

The  computer  code  results  for  the  NPS  HRGK  design  case  are  shown  in 
Figure  23.  The  results  are  for  the  maximum  launch  range  case  (6  km)  as  this  is  the  most 
stressing  case.  All  other  input  parameters  are  set  as  in  previous  analyses.  The  target  speed 
is  uniformly  distributed  between  20  to  60  miles  per  hour.  The  results  show  that  a  95% 
probability  of  having  the  target  in  the  field  of  view  at  time,  can  be  obtained  with  fields  of 
view  slightly  smaller  than  those  that  would  be  required  for  the  same  results  against 
stationary  targets.  For  example,  with  a  seeker  acquisition  range  of  1.5  km,  a  ±25°  seeker 
would  provide  a  95%  probability  of  acquisition  against  the  probabilistic  target  model.  But, 
against  a  stationary  target,  a  +29°  seeker  would  be  required  for  the  same  confidence  (Figure 
17,  page  40).  The  probabilistic  moving  target  model  results  in  a  smaller  field  of  view 
because  of  the  possibility  of  the  target  moving  in  the  same  direction  as  the  rocket  cross¬ 
track  error  is  allowed.  This  possibility  increases  the  probability  of  having  the  target  within 
the  field  of  view,  even  with  a  slightly  smaller  field  view. 

D.  COMMENTS  AND  OBSERVATIONS 

In  concluding  this  chapter  on  preliminaiy  sizing  analyses,  several  observations  are 
appropriate.  These  include  (1)  a  discussion  of  the  benefits  and  data  that  came  from  the 
computer  models,  (2)  a  word  of  caution  about  the  limitations  of  the  models,  and  (3)  a 
listing  of  additional  work  that  could  be  done  to  expand  the  usefulness  of  the  preliminary 
sizing  models. 


'  Actually  the  probability  of  the  target  being  within  the  cross-track  limits  of  the  seeker’s  field  of  view. 
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1.  Additional  Data  from  Modeling 

The  primary  purpose  of  the  computer  models  discussed  in  this  chapter  was  to  give 
the  designer  a  better  understanding  of  how  the  system  requirements  affect  different  design 
specifications.  However,  the  codes  not  only  facilitate  the  “flow-down”  of  requirements,  but 
they  also  provide  additional  insights  into  the  design  of  the  system.  In  the  case  of  the  NPS 
HRGK  design,  several  observations  were  made  from  the  use  of  the  preliminary  sizing 
codes.  The  following  paragraphs  site  some  examples. 

From  the  length,  weight,  and  center  of  gravity  analysis  it  was  determined  early  in 
the  design  that  the  launcher  center  of  gravity  constraint  is  not  binding  for  any  design  within 
the  maximum  weight  and  length  constraints.  Thus,  the  kit’s  center  of  gravity  can  be  placed 
to  provide  the  best  combination  of  performance  and  packaging  without  concern  for  the 
launcher’s  center  of  gravity.  In  a  general  case,  a  designer  can  use  this  tjqie  of  knowledge  to 
move  through  design  iterations  without  making  calculations  regarding  the  non-binding 
constraints.  Then,  as  the  final  design  is  approached,  calculations  can  be  made  to  verify  the 
preliminary  assessments. 

Another  good  example  of  additional  information  gained  from  the  models  is  the 
minimum  required  seeker  acquisition  range.  The  required  seeker  range  was  inadvertently 
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discovered  by  trying  several  seeker  acquisition  ranges  in  the  maneuverability  analyses.  The 
final  seeker  range  specification  could  be  derived  based  on  any  of  the  three  target  motion 
types  (stationary,  worst-case,  or  uniformly  distributed). 

The  maneuverability  and  field  of  view  analyses  provided  other  general  insights  into 
the  effect  of  the  system  requirements.  For  example,  from  the  analyses  it  became  apparent 
that  the  minimum  launch  range  requirement  defines  the  maneuverability  specification  and 
that  the  maximum  launch  range  requirement  drives  the  seeker  field  of  view  specification. 
This  insight  is  useful  in  determining  design  tradeoffs  with  relaxed  system  requirements. 
For  example,  if  the  minimum  range  requirement  for  the  NFS  HRGK  were  increased  to  2 
km,  the  effect  on  the  seeker  would  be  minimal,  but  the  affect  on  the  thruster  design  would 
be  significant.  Knowing  this  the  designer  can  determine  the  value  of  relaxing  the  system 
requirement. 

2.  Limitations  of  Models 

Despite  the  usefulness  of  the  codes  discussed  in  this  chapter,  some  caution  must  be 
applied  when  using  these  and  other  mathematical  or  computer  models.  Model  outputs  are 
sometimes  skewed  or  erroneous  because  of  (1)  the  nature  of  the  model  with  its  simplifying 
assumptions  or  (2)  because  of  the  model’s  sensitivity  to  incorrect  inputs.  Examples  of  both 
these  error  sources  are  abundant  in  the  NFS  HRGK  analyses. 

For  example,  both  the  maneuverability  and  the  field  of  view  analyses  used  a  two- 
dimensional  model  for  a  three  dimensional  problem.  Implied  in  these  models  is  the 
assumption  that  the  azimuth-related  maneuverability  and  field  of  view  requirements 
outweigh  the  elevation-related  requirements.  In  reality,  the  HRGK  must  execute  significant 
maneuvers  in  both  yaw  and  pitch  (at  the  same  time)  to  prosecute  a  target  as  it  appears  in  the 
seeker’s  field  of  view.  Furthermore,  the  field  of  view  sizing  may  well  be  defined  by  the 
requirement  to  look  down  at  a  target  rather  than  the  requirement  to  look  to  the  left  or  right 
for  a  target.  These  possibilities  are  not  accounted  for  in  the  two-dimensional  model. 

Another  example  of  potentially  incorrect  modeling  concerns  the  appropriateness  of 
the  target  motion  model.  The  constant  heading  and  speed  with  uniform  distributions  may  be 
an  appropriate  target  motion  model  for  short  time  of  flight  attacks  against  sluggish  ground 
targets.  However,  this  type  of  target  motion  model  would  not  be  appropriate  in  a  surface- 
to-air  scenario  where  the  target  may  have  evasive  maneuver  capability.  In  general,  any 
target  motion  model  must  allow  for  a  reasonable,  worst-case  evasive  maneuver. 
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Figure  24.  Sensitivity  Analysis  of  the  Steady-Turn  Maneuver  Model  for  Thruster  Sizing, 
(a)  Cm  ±  100%  and  (b)  Xcp  ±  3%. 

Finally,  Figure  24  illustrates  the  importance  of  sensitivity  analysis  and  the  potential 
for  erroneous  results  due  to  incorrect  inputs.  The  figure  shows  data  from  the  control 

thruster  sizing  analysis  with  variations  in  the  steady- turn  model  parameters.  The  plot  on  the 
left  (a)  shows  the  effect  of  ±100%  changes  in  Cw  .  Doubling  or  zeroing  the  parameter 

produces  only  minor  changes  in  required  control  thrust.  Clearly,  the  thrust  required  is  very 
insensitive  to  changes  in  On  the  other  hand,  the  plot  on  the  right  shows  the  impact  of 

±3%  changes  in  Xcp.  The  small  changes  in  this  input  parameter  affect  the  required  control 
thrust  significantly. 

Another  point  of  interest  illustrated  in  the  above  example  is  the  importance  of 
parameter  selection.  In  the  NFS  HRGK  case,  the  parameter  X(.p  is  measured  from  the  nose 
of  the  guided  rocket  and  varies  between  28.5  and  31.8  inches  as  a  function  of  Mach 
number.  Small  percentage  changes  in  Xcp  result  in  large  percentage  changes  in  the  rocket’s 
static  margin  (the  small  distance  between  the  x^p  and  the  x^g).  In  fact,  the  ±3%  change  in 
Xcp  produces  between  a  ±22.7%  and  ±86.5%  change  in  static  margin.  Although  either  X(.p 
or  static  margin  could  be  used  as  model  inputs,  the  model  sensitivity  to  percent  changes  in 
the  two  parameters  would  vary  significantly. 

3.  Additional  Work 

The  computer  models  developed  by  the  author  could  be  modified  to  produce  new  or 
refined  models.  The  addition  of  probabilistic  target  motion  to  the  turn  rate  and  turn  radius 
code  would  provide  useful  data.  The  analysis  of  elevation  plane  considerations  would  also 
be  helpful  in  the  guided  rocket  design.  A  model  that  could  estimate  both  the  required 
vertical  field  of  view  and  pitch  axis  maneuverability  would  help  in  verifying  the  thruster 
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IV.  SIX-DEGREE  OF  FREEDOM  SIMULATION 


This  chapter  addresses  the  use  of  detailed  computer  modehng  in  the  missile  design 
process.  The  chapter  describes  a  six-degree  of  freedom  (6DOF)  simulation  developed  by 
the  author  for  the  analysis  of  the  NFS  HRGK.  The  HRGK  6IX)F  is  a  major  modification 
to  a  simulation  developed  by  Professor  Robert  (Gary)  Hutchins,  Electrical  and  Computer 
Engineering  Department,  Naval  Postgraduate  School,  and  his  1996  “Navigation,  Missile, 
and  Avionics  Systems”  class.  The  simulation  was  developed  in  Simulink  using  Mailab 
functions  to  perform  the  simulation’s  computations.  The  listings  of  the  Matlab  functions 
are  included  in  Appendix  C. 

Currently,  the  6DOF  code  requires  more  than  20  minutes  to  simulate  a  15  second 
rocket  flight^  To  speed  up  the  simulation,  some  of  the  simulation  features  are  not  yet  fully 
implemented,  but  the  author  is  in  the  process  of  converting  the  simulation  to  a  more 
efficient  computation  structure.  This  structure  uses  the  native  Simulink  “S-function”  and 
was  introduced  to  the  author  by  Mr.  Allen  Robins;  Dynamics  and  Controls  Section,  Naval 
Air  Warfare  Center  Weapons  Division,  China  Lake,  California.  (Appendix  B  briefly 
describes  the  different  computation  structures  the  can  be  used  in  Simulink.)  Preliminary 
tests  by  the  author  have  shown  that  the  use  of  the  S-function  structure  can  cut  simulation 
run  times  by  a  factor  of  nearly  30. 

This  chapter  begins  with  an  overview  of  the  simulation  architecture  and  then 
addresses  the  functions  of  the  major  simulation  components.  A  section  of  the  chapter  is 
devoted  to  the  development  of  the  guided  rocket  equations  of  motion.  And  another, 
presents  simulation  results  used  in  the  design  of  the  NPS  HRGK.  The  chapter  ends  with  a 
few  observations  concerning  the  use  of  detailed  simulations  in  the  missile  design  process. 

A.  SIMULATION  OVERVIEW 

The  HRGK  6DOF  simulation  is  designed  to  be  used  with  an  initialization  script  file. 
This  Matlab  file  is  executed  prior  to  running  the  simulation.  It  sets  up  simulation 
parameters  such  as  vectors  with  the  rocket’s  aerodynamic  coefficients  (functions  of  Mach) 
and  time  profiles  for  mass  properties  and  the  main  rocket  thrust  The  file  also  sets  the  initial 
conditions  for  all  the  simulation  states.  When  run,  the  simulation  uses  a  user-selectable 

'  Silicon  Graphics,  Inc.  Indigo  II  workstation 
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integration  scheme  to  simulate  the  rocket  flight.  The  simulation  writes  several  output 
vectors  to  the  Matlab  workspace.  These  outputs  include  the  flight  time,  target  miss 
distance,  the  guided  rocket’s  states  (f/,  V,  W,  P,  Q,  R,  North,  East,  Altitude,  (f),  6,  and 
Iff),  and  their  time  derivatives  with  respect  to  the  appropriate  coordinate  frames. 

Figure  25  shows  a  top-level  block  diagram  of  the  HRGK  6DOF.  Figure  26  and 
Figure  27  respectively  show  detailed  views  of  the  “Rocket  Dynamics”  and  the  “Parameters 
and  Coefficients”  blocks  as  they  are  currently  implemented. 


Figure  25.  HRGK  6DOF  Simulation  Overall  Block  Diagram.  (Dashed  blocks  are  not  yet  fully 
implemented.) 


The  blocks  in  Figure  25  with  dashed  lines  are  not  yet  fully  implemented.  For 
example,  the  thruster  control  inputs  are  not  currently  generated  based  on  seeker  inputs,  but 
rather,  the  thruster  commands  must  be  “hard  coded”  into  the  controller  block. 

Inside  the  Rocket  Dynamics  block,  the  time  derivatives  of  the  rocket’s  states  are 
computed  by  a  Matlab  function.  The  state  derivatives  are  then  integrated  and  feedback  into 
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the  function  and  other  parts  of  the  simulation.  The  rocket’s  aerodynamic  and  mass  property 
parameters  are  determined  using  lookup  tables  in  the  Parameters  and  Coefficients  block. 


Figure  26.  Detailed  View  of  the  Rocket  Dynamics  Block. 
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Figure  28  shows  the  overall  architecture  for  the  simulation  currently  under 
development.  Most  of  the  components  of  this  architecture  are  described  in  the  following 
section. 
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Figure  28.  Future  6DOF  Architecture. 


B.  SIMULATION  COMPONENTS 

The  simulation  shown  schematically  in  Figure  28  is  made  up  several  components. 
These  components  are  very  briefly  described  in  the  following  subsections. 

1.  Environment 

The  environment  block  simulates  the  environment  in  which  the  rocket  operates. 
This  block  produces  the  simulation  time  and  provides  constant  gravity.  It  also  computes  the 
atmospheric  properties  needed  to  calculate  the  Mach  number  and  dynamic  pressure.  The 
block  uses  the  Dryden  wind  gust  model  to  produce  wind  gusts  in  pitch  and  yaw. 

2.  Missile 

The  missile  block  simulates  the  dynamics  of  the  guided  rocket.  The  block  has 
several  sub-components.  These  blocks  are  divided  to  align  as  closely  as  possible  with 
actual  subsystems.  The  missile  block  uses  external  signals  and  data  to  produce  realistic 
seeker  and  IMU  signals.  From  these  signals  control  inputs  are  determined  and  applied  to 
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the  forces  and  moments  sub-component.  The  forces  and  moments  component  computes  all 
the  body  frame  forces  on  the  rocket.  From  these  forces  and  moments  the  dynamics  of  the 
body  can  be  determined  by  the  flat-earth  equations  of  motion. 

3.  Target  Dynamics 

The  target  block  provides  the  target  location  and  motion  to  the  relative  geometry 
block.  This  data  allows  the  line  of  sight  (between  the  rocket  and  the  target)  and  the  line  of 
sight  rate  data  to  be  calculated. 

4.  Relative  and  Absolute  Geometry 

The  relative  geometry  block  uses  the  target  and  rocket  dynamics  to  solve  the 
navigation  equations  (described  in  the  following  section).  From  the  navigation  equations, 
the  absolute  and  relative  position  of  the  rocket  and  target  can  be  can  be  computed.  The  line 
of  sight  and  line  of  sight  rates  are  also  provided  by  the  geometry  block. 

C.  EQUATIONS  OF  MOTION 

This  section  contains  a  brief  derivation  of  the  flat-earth  equations  of  motion  as  they 
apply  specifically  to  the  HRGK  6DOF.  The  derivation  is  adapted  from  Stevens  and  Lewis 
[Ref.  15]. 

1.  General  Matrix-Form  Equations 

Starting  with  generic  flat  earth  equations  of  motion: 

F 

Vg  =  -QjV  g  +  Bggp  -I-  -S-  (Force  Equation) 

m 

(Moment  Equation) 
O  =  E(<I>)©  Q  (Attitude  Equation) 

J?NED  -  ^B^ B  (Navigation  Equation) 

where  the  is  the  velocity  vector,  [U  V  W]^;  go  is  the  gravity  vector  at  the  surface  of  the 

earth  (adjusted  for  the  earth  rotation),  [0  0  9.81m/sf ;  F^  is  the  force  vector,  [F^  F^  FJ’'; 
(Ob  is  the  angular  rate  vector,  [P  Q  R]^;  Tg  is  the  torque  vector,  [L  M  N]"^;  Hg,  to  be 
defined  later,  is  a  vector  to  account  for  the  gyroscopic  effects  of  the  spinning  rocket  motor; 
^  is  the  Euler  angle  vector,  0  is  the  north-east-down  frame  position  vector, 

tPw  Pe  PdI  ’  ^  is  the  body  frame  angular  rates  cross-product  matrix. 
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Qg  = 


0  -R  Q' 
R  0  -P 
-Q  P  0 


is  the  rotation  matrix  from  the  North-East-Down  frame  to  the  body  frame, 


'1 

0 

0 

COS0 

0 

— sin0 

cosy/ 

sin^ 

O' 

Bg  = 

0 

COS0 

sin^ 

0 

1 

0 

-siny/ 

cosy/ 

0 

0 

— sin^ 

COS0 

sin0 

0 

COS0 

0 

0 

1 

J  is  the  inertia  matrix  (with  symmetry  assumed  about  the  x-axis), 


7  = 


Jx  0 
0  7v 


0 

0 


0  0  7. 


and  the  Euler  angle  function  is 


£(®)  = 


1  tan0sin0  tan^cos^ 

0  cos<l>  —sm<p 

Q  sin^  cos^ 


COS0  COS0 

The  state  vector  is  Pned^^-  Now  the  forces  and  moments  need  to  be  defined 

for  the  HRGK  case. 

2.  Forces  and  Moments 

The  aerodynamic  forces  and  moments  are  functions  of  Mach  number  and  can  be 
expressed  by  the  coefficients  and  their  derivatives.  The  following  summarizes  the  forces 
and  moments  computed  in  HRGK  6DOF: 


f  2  \ 

F^~T-U 

(Axial  Force) 

Fy  =  Ty+ 

(Side  Force) 

Pz  ~^z~ 

(Down  Force) 

(C,  d  \ 

(Rolling  Moment  on  HRGK) 
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(Pitching  Moment) 


M  =  T^XT  +  qAd 


N  =  TyXT  +  qAd 


CN^OLXcp-^-TjfQ 

ZU  J 


(Yawing  Moment) 


where  T  is  the  rocket  motor  thrust;  and  are  they  and  z  components  of  the  reaction 
thruster  forces,  =  r,cos0  +  r2Cos(^2/37i)  +  T^cos(<lH-4/3n)  and  =  r,sin^  + 
r2sin(^2/37c)  +  r3sin(^f4/37r);  and  is  the  rolling  torque  exerted  on  the  HRGK  by  the 
spinning  rocket  motor, 

Tr  = 

where  fj.  is  the  coefficient  of  friction  for  the  bearings,  is  the  axial  moment  of  inertia  for 
the  rocket  alone,  and  Pj^  is  the  rocket  motor’s  spin  axis  angular  acceleration  (adapted  from 
Ref.  [11]).  The  state  has  to  be  added  to  the  state  vector,  and  the  gyroscopic  affect  of  the 
spinning  rocket  is  added  to  the  moment  equation.  The  effect  can  be  accounted  for  by  adding 
the  vector  to  the  right-hand  side  of  the  moment  equation,  where  H,^  is 

MT 

J  J 

y  y 


3.  HRGK  Scalar  Flat-Earth  Equations 

From  the  general  matrix  equations,  the  HRGK-specific  scalar  form  of  the  equations 
can  be  written  in  terms  of  the  defined  forces  and  moments. 

a.  Force  Equations 

U  =  RV-QW-g^sme  +  ^ 

m 

V  =  -RU  +  PW  +  g^sirKpeosO  +  — 

m 

W  =  QU  -  PV  +  gQCOs<l>cosd 

m 


b.  Moment  Equations 
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Q  = 


1-^ 


M  4 

PR  + - 

J  J 

y  y 


R  = 


^  J  ^ 

ii-l 

\^y  J 


N  J 

PQ  +  y  +  ^QP, 


c.  Kinematics  Equations 

^  =  P  +  tan0(<2sin0  +  /?cos0) 

9  =  Qcos(l>-  Rsin<j> 
,cos<}> 


\jf  =  Qsin^  +  R 


cosd 


d.  Navigation  Equations 

Pf^  =  Ucosdcosy/  +  V{sm<l>sindcosy/ -cos(}>smY)+W{cos(l>sindcosY  +  sin<l)siny/) 
Pfj  =  Ucos0sinY+V{sin<psindsin\if  +  cos<t>cosy/)  +  W{cos(t>&ir\6smy/-sin(j)cosiit) 
Pq  = -Usin9  +  Vsin^cos9  +  Wcos(l>cosd 


D.  NPS  HRGK  RESULTS 

The  force  and  moment  equations,  as  well  as  the  equations  of  motion,  are 
implemented  in  a  Matlab  function  which  is  used  by  the  HRGK  6DOF  simulation.  The 
following  subsections  show  two  examples  of  how  the  simulation  was  used  in  the  design. 

1.  Ballistic  Flight  and  Maximum  Range  (3DOF) 

Early  in  during  the  NPS  HRGK  design  process,  the  HRGK  simulation  was 
simplified  to  only  include  the  three  dimensions  of  the  rocket’s  longitudinal  dynamics.  (The 
gyroscopic  effects  of  the  spinning  rocket  were  also  ignored.)  Using  the  simplified  3DOF, 
the  results  shown  in  Figure  29  were  obtained.  These  results  confirmed  that  the  HRGK- 
equipped  rocket  could  meet  the  time  to  target  and  the  maximum  range  requirements. 

2.  Pulse  Width  Modulated  Turn  Rate  (6DOF) 

Later  in  the  design,  the  affect  of  thruster  pulse  width  modulation  was  investigated 
using  the  HRGK  6DOF.  Two  simulations  were  run.  The  first  simulated  the  missile  flight 
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Figure  29.  HRGK-Equipped  Rocket  Ballistic  Trajectories.  Ground  launch  at  sea 
level  with  various  launch  elevation  angles.  Solid  lines  are  for  cold  rocket  motor, 
dashed  lines  are  for  hot  rocket  motor,  and  flight  times  are  indicated. 

with  a  constant  thruster  force  from  the  appropriate  side  only.  The  second  used  a  pulsed 
sequence  of  all  three  thrusters,  but  with  a  longer  pulse  on  the  appropriate  side.  The  results 
for  the  two  runs  are  shown  in  Figure  30.  From  the  analysis,  it  was  determined  that  constant 
thrust  does  not  generate  satisfactory  maneuvers.  (In  fact,  in  the  simulation,  the  rocket 
began  “cartwheels”  without  making  substantial  changes  in  cross-track).  On  the  other  hand, 
the  pulsed  control  generated  an  adequate  maneuver  while  maintaining  a  somewhat 
reasonable  side  slip  angle,  p. 


Figure  30.  Sustained  and  Pulsed  Thruster  Maneuvers.  Flight  paths  and  headings  indicated  for 
a  0.8  second  maneuver. 
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E.  OBSERVATIONS 


As  designs  become  more  refined  and  leave  the  conceptual  phases,  high  fidelity 
simulations  becomes  more  important.  In  the  case  of  the  NFS  HRGK,  the  6DOF  simulation 
was  helpful  in  performing  some  detailed  analyses.  Nevertheless,  the  time  required  per 
simulation  run  precluded  the  use  of  the  simulation  in  many  cases. 

Fast,  modular  simulations  that  can  be  easily  modified  for  various  designs  are 
desirable  design  tools.  In  many  cases,  lower-order  simulations  (3DOF  or  5DOF)  meet  all 
the  needs  of  the  designer.  For  example,  the  3DOF  version  of  the  HRGK  simulation  had 
sufficient  fidelity  to  be  used  in  the  maximum  range  and  time  to  target  analyses. 
Unfortunately,  other  analyses  required  a  6DOF  so  that  the  roll  of  the  seeker  and  the 
gyroscopic  effects  of  the  spinning  rocket  could  be  observed. 

If  the  restmcturing  of  the  HRGK  simulation  provides  the  expected  improvement  in 
run  time,  this  simulation  would  become  a  much  more  useful  design  tool.  With  run  times  of 
only  a  minute  or  so,  several  design  iterations  could  be  evaluated  quickly. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  purpose  of  this  thesis  is  to  provide  documentation  for  the  several  computer 
codes  developed  by  the  author  during  the  design  and  analysis  of  the  NFS  HRGK.  The 
thesis  is  also  intended  to  provide  insights  concerning  the  use  of  modeling  in  the  design  and 
analysis  processes.  Several  germane  points  can  be  made  in  conclusion. 

First,  models  with  differing  levels  of  complexity  are  appropriate  during  different 
phases  of  the  design.  Models  used  in  the  analysis  of  conceptual  and  preliminary  designs  are 
usually  the  simplest  while  those  models  used  in  design  refinement  are  more  complicated. 
COEA-related  modeling  is  essential  throughout  the  design  process  to  ensure  that  the  design 
meets  the  system’s  cost  effectiveness  goals. 

Second,  simplistic  models  can  be  used  to  ensure  the  flow  down  of  system 
requirements  into  design  specifications.  In  this  way,  the  models  help  ensure  a  design 
driven  by  the  requirements.  Often  these  simple  models  provide  great  insight  into  additional 
aspects  of  the  design  problem  allowing  several  critical  specifications  to  be  defined  in  the 
earliest  phases  of  the  design. 

Third,  modeling  in  the  design  process  should  include  the  appropriate  mix  of 
probabilistic  and  worst-case  conditions.  These  two  types  of  models  can  provide  very 
different  results  which  must  be  judiciously  evaluated. 

Finally,  knowing  the  sensitivity  of  a  model  to  variations  in  parameters  is  critical  to 
the  proper  use  of  the  model  in  decision  making.  Sensitivity  analyses  or  “what-if  ’  studies 
must  be  performed  to  ensure  the  modeling  is  providing  reliable  data. 

In  addition  to  these  concluding  points  several  recommendations  for  foUow-on  work 
are  also  appropriate.  First,  the  design  and  analysis  of  the  NFS  HRGK  would  have  been 
facilitated  by  simple  Matlab  codes  that  allow  the  early  aerodynamic  sizing  of  the  kit.  Such 
a  program  would  be  helpful  in  evaluating  the  effect  of  nose  shape  and  canard  size  and 
shape  on  the  range  of  the  HRGK-equipped  rocket 

Additional  work  is  also  needed  to  include  a  probabilistic  target  motion  model  in  the 
maneuverability  analysis  code.  The  development  of  models  that  evaluate  the  required  pitch 
maneuverability  and  the  required  elevation  field  of  view  would  also  be  useful. 

Finally,  the  modification  of  the  HRGK  simulation  to  make  it  a  generic  bDOF 
simulation  would  be  very  useful  for  future  missile  design  work.  The  Simulink  architecture 
can  offer  fast  modular  simulations  that  can  be  easily  modified  to  meet  future  users’  needs. 
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APPENDIX  A.  MISSILE  AERODYNAMICS  NOMENCLATURE 


Many  of  the  aerodynamic  coefficients,  stability  derivatives,  and  various  other  terms  and 
symbols  used  in  this  thesis  are  defined  below.  The  nomenclature  closely  follows  that  used 
by  Chin  [Ref  16]  with  some  adaptations  like  that  used  by  Blakelock  [Ref  17]  and  by 
Howard  [Ref  18].  It  is  worth  noting  two  significant  differences  between  terms  associated 
with  missiles  and  those  associated  with  traditional  aircraft.  First,  most  missile  coefficients 
are  referenced  to  the  missile's  diameter  or  its  cross-sectional  area  (normal  to  the 
longitudinal  axis)  rather  than  a  wing  chord  or  a  planform  area.  Second,  normal  and  axial 
force  coefficients,  and  C^,  (aligned  to  the  body-axis  system)  are  used  rather  than  the 
lift  and  drag  force  coefficients.  Cl  and  Cp,  (in  the  stability  axis  system).  When  the 
missile's  angle  of  attack  (a)  is  small  there  is  very  little  difference  between  the  two  systems, 
but  for  significant  angles  of  attack,  the  relationship  given  in  the  following  equations  should 
be  used. 

Ca  =  Cj)  cos  a  -  Cl  sin  a 
Cn  =  Cl  cos  a  -  Cj)  sin  a 

Reference  Dimensions 

A  or  Cross-sectional  area  of  missile  body  (5.94  in^) 
d  Diameter  of  missile  body  (2.75  in) 

L  Missile  length  (68.75  in) 

Xgg  Distance  from  nose  tip  to  center  of  gravity  on  missile  (27  to  3 1 .7  in) 

Flight  Conditions 

a  angle  of  attack 

P  side  slip  angle 

M  Freestream  Mach  number 

V,  Freestream  air  speed 

p  Freestream  air  density 

q  Dynamic  pressure  ( j  pV^ ) 
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Aerodynamic  Coefficients* 


Cl 

Cn 

Cm 

Ca 

Cd 


Coefficient  of  lift, 


Lift 


Normal  force  (N)  coefficient, 


N 

qSn 


Pitching  moment  (M)  coefficient. 


Axial  force  (A)  coefficient. 


M 

qSnd 


Drag  force  (D)  coefficient. 


D 

q^n 


Cdq  Drag  force  coefficient  at  zero  lift 

Cdj^  Coefficient  of  drag  due  to  normal  force  (D^), 

Static  and  Dynamic  Derivatives* 


Pn 

qS;r 


dC. 


Cn  Normal  force  coefficient  per  radian  of  angle  of  attack,  — ^ 
“  da 

Cm„  Pitching  moment  coefficient  per  radian  of  angle  of  attack. 


dC 


M 


da 

C^^  Pitching  moment  coefficient  per  radian/sec  of  pitch  rate,  f  — 
Missile  “States” 


U,  V,  W  Body  axis  velocities  (x,  y,  and  z-axes) 

P,  Q,R  Body-axis  angular  rates  (about  the  x,  y,  and  z-axes) 

(t>,  0,  i/r  Euler  Angles  (roll,  pitch,  yaw) 

PfpPpPD  Earth  Positions  (North,  East,  Down) 
h  height  (-P^) 


*  Directional  coefficients  and  derivatives  (e.g.-,  Cy,  C;;^,  Cy^,  C  Cj;j  )  are  assumed  equal  to  their 
longitudinal  counterparts  (e.g.;  C^,,  Cm,  Cn„,  Cwa,  Cm,;  respectively). 
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APPENDIX  B.  NOTES  ON  MATLAB  AND  SIMULINK 


A.  MATLAB 

Matlab  is  a  high-performance  numeric  computation  and  visualization  computer 
software  package.  It  is  essentially  a  high-level  programming  language  that  is  especially 
tailored  for  the  manipulation  of  matrices,  (see  Matlab  Reference  Manual  and  Matlab 
User’s  Guide) 

Matlab  commands  can  either  be  invoked  in  an  interactive  mode  at  the  keyboard  or 
through  Matlab  script  files.  Another  type  of  file  is  the  Matlab  function.  Both  the  script 
and  function  files  are  labeled  with  an  “m”  extension  at  the  end  of  the  file  name  and  are 
commonly  called  M-files. 

A  script  M-file  can  run  as  a  stand  alone  “batch”  file  invoking  Matlab  commands 
when  the  file  is  executed.  On  the  other  hand,  a  function  M-file  is  called  from  within  script 
files  or  from  the  keyboard.  A  function  can  be  called  with  or  without  arguments  and  can 
return  computed  values  when  completed. 

Variables  declared  or  set  in  a  script  file  or  at  the  keyboard  are  available  throughout 
the  Matlab  “workspace.”  These  variables  are  common  and  can  be  accessed  by  other  script 
files.  On  the  other  hand,  variables  declared  in  a  function  M-file  are  local  variables  and  are 
not  available  in  the  workspace.  Furthermore,  functions  do  not  have  access  to  workspace 
variables  unless  the  variables  are  passed  to  the  function  or  are  declared  as  “global”  in  both 
the  workspace  and  in  the  function. 

The  standard  syntax  used  in  the  first  line  of  a  function  M-file  is 

function  y  =  funciu) 

where  “function”  is  the  Matlab  function  indicator,  y  is  the  variable  that  wiU  be  passed 
back,  u  is  the  input,  and  func  is  the  name  of  the  function  (func.m). 

Files  that  are  written  in  C  or  FORTRAN  can  be  converted  into  executable  Matlab 
files  with  a  “.mex”  extension  by  using  the  “cmex”  or  “fmex”  commands,  (see  Matlab 
External  Interface  Guide)  Mex-files  have  the  advantage  of  running  faster  than  normal  M- 
files. 
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B.  SIMULINK 


Simulink  is  a  dynamic  system  simulation  software  package  and  an  extension  to 
Matlab.  Simulink  offers  many  ready-made  blocks  that  can  be  assembled  graphically  to 
model  a  dynamic  system.  Two  special  Simulink  blocks  allow  tremendous  flexibility  in 
modeling  systems,  (see  Simulink  User’s  Guide) 

The  first  is  the  Matlab  function  block.  This  block  when  placed  in  a  block  diagram 
will  execute  the  commands  in  the  indicated  M-file  function.  So  for  example,  the  equations 
of  motion  for  a  system  could  be  placed  in  a  M-file  function.  This  function  would  have 
inputs  from  simulink  that  include  the  states  and  the  control  input.  The  function  would  then 
output  the  states’  time  derivatives  which  could  then  be  integrated  and  then  passed  back  the 
Matlab  function. 

This  method  runs  extremely  slow! 

The  second  more  useful  block  is  the  S-function  block.  This  block  calls  either  a  M- 
file  or  a  Mex-file  that  is  specifically  stmctured  to  run  efficiently  in  Simulink  (see  Simulink 
1.3  Release  Notes).  Chapter  4  of  Simulink  1.3  Release  Notes  outline  the  procedure  for 
creating  Mex-files  from  C  codes  for  use  in  fast  mnning  Simulink  models. 

Additionally,  simulink  models  that  do  not  contain  M-file  references  (that  use  Mex- 
files  in  the  S-function  blocks)  can  be  turned  into  stand-alone,  accelerated,  executable  code 
using  Simulink’s  Real-Time  Workshop  (see  Real-Time  Workshop  User’s  Guide). 
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APPENDIX  C.  MATLAB  CODE  LISTINGS 


This  appendix  provides  listings  of  the  following  Matlab  computer  codes  which 
were  discussed  in  this  thesis: 

•  Single-Shot  Probability  of  Kill 

•  Length,  Weight,  and  Center  of  Gravity 

•  Turn  Rate  and  Radius 

•  Steady-Tum  Maneuver 

•  Seeker  Field  of  V iew 

•  6DOF  Codes 
-Initialization  Code 
-Aerodynamics  Data 
-Equations  of  Motion 
-Line  of  Sight 
-Mach 

The  user-defined  inputs  in  the  codes  are  denoted  with  bold-face  type. 


C.  SINGLE-SHOT  PROBABILITY  OF  KILL 

%  Single  Shot  Pk  using  "Cookie  Cutter"  methodology 

%  W.  Mark  Wonnacott,  Naval  Postgraduate  School,  May  1997,  corrected/revised  Aug  1997 

%  Terminal  Dive  Angle  (deg  relative  to  horizontal) 

dive  «  20; 

%  CEPs  to  be  used  in  computations  (for  a  single  CEP  set  CEP_min  equal  to  CEP_max) 
CEP_min  ■  .5;  %  minimum  CEP  (m)  to  be  used 

CEP_max  s  loO;  %  maximum  CEP  (m)  to  be  used 

%  Target  Size  (m) ,  Prob.  of  Kill  given  Hit,  &  Desired  Prob.  Kill 

L  =  5.5;  H  =  2.0;  W  =  2.5;  Pkh  =.90;  PkD  =.95;  trgt= '  Light  Arm.  (AirDef )  •  ; 

L  =  5.5;  H  =  2.0;  W  =  2.5;  Pkh  =.80;  PkD  =.60;  trgt=' Light  Armor  (APC)  ’; 

L  =  8.0;  H  =  2.5;  W  =  2.5;  Pkh  =  1;  PkD  =.50;  trgt=* Large  Truck'; 

L  =  14.;  H  =  1.5;  W  =  4.5;  Pkh  =,70;  PkD  =.75;  trgt='Patrol  Boat '  ; 

%  Elevation  and  Crossrange  Projected  Dimensions  (l:Head  on  &  2:Broadside) 

Yl  =  H*cos (dive*pi/180) +L*sin (dive*pi/180) ; 

XI  =  W; 

Y2  =  H*cos (dive*pi/180) +W*sin (dive*pi/180) ; 

X2  =  L; 


%CEPs  in  desired  range 

if  CEP_min<=0,  CEP_min=eps;  end  %CEP_min  set  to  allowable  small  number 
if  CEP_max<=0,  CEP_max=eps;  end  %CEP_max  set  to  allowable  small  number 
if  CEP_min>CEP_max,  temp=CEP„min;  CEP_max=CEP_min;  CEP_min=temp;  end 
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CEPs  =  logspace (logic  (CEP_jnin) , loglO  (CEP_mx) )  ; 

sigmas  =  CEPs/1.1774;  %corresponding  sigmas,  1-D  miss  distances  (m) 

S  =  sigmas*sqrt  (2)  *2;  %norTnalizd  half  sigmas  for  Matlab 

%  Probability  of  Hit  (equal  (50%)  prob.  of  broadside  or  head-on  encounter) 

Ph  =  0.5*(erf  (Yl./S)  .*erf  (Xl./S) )  +  0.5*  (erf  (Y2  . /S)  .  *erf  (X2  . /S) )  ; 

%  Probability  of  Single  Shot  Kill 
PkSS  =  Pkh*Ph; 

%  Number  of  Weapons  Needed  to  Achieve  Desired  Prob.  of  Kill 
N  =  log(l-PkD) ./log(l-PkSS); 

figure  (1)  %  Single-Shot  Probability  of  Kill 

semilogx  (CEPs, PkSS, '-')  ,  grid  on 
axis  ( [CEP_inin,  CEP_inax ,0,1]) 

xlabel  ( '  CEP  (m)  ' ) ,  ylabel  ( '  Single  Shot  Probability  of  Kill  ’ ) 

figure  (2)  %  Number  of  Weapons  per  Kill 

loglog  (CEPs,N, '-' ) ,  grid  on 
axis  ( [CEP_min, CEP__max,  1 , 1000] ) 

xlabel ( 'CEP  (m)  ’ )  /  ylabel ( [  'Number  of  Rockets  for  Pk  =  '  ,num2str  (PkD)  ] ) 

figure  (3)  %  Relative  Cost  for  Equal  Cost /Kill 

Nc  =  max(l,ceil (N) )  ; 

[X,  Y]  =stairs  (CEPs , Nc  ( 1 )  .  /Nc)  ; 
plot  (X,Y,  ’-•) ,  grid  on 
axis  ( [CEPjmin,  CEP_inax,  0,1]) 

xlabel  ('CEP  (m)  ' ) ,  ylabel  ( 'Relative  Weapon  Cost' ) 


D.  LENGTH,  WEIGHT,  AND  CENTER  OF  GRAVITY 

%  GUIDED  Rocket  Center  of  Gravity  as  a  Function  of  Kit  Weight 
%=========:==========================:========================= 

%  W.  Mark  Wonnacott,  Naval  Postgraduate  School,  May  1997 


%%%%  Design  Space  Paramters 

Wkl  a  linspace ( 1, 8} ;  % 

Lkl  a  [5,10,13.625,Lmax-Lr] ;  % 

Dkl  =  [.05, .07, .0865, .1, .IS] ;  % 

Cgkla  1/3;  % 

dla  a  2.75;  % 

areaa  pi* (dia/2 ) ^2;  % 


as  summed  kit  mass  range  (Ibm) 

specific  kit  lengths  (in) 

specific  kit  densities  (lbm/in^3) 

kit  Cg  from  nose  (fraction  of  kit  length) 

kit  diameter  (in) 

kit  cross  section  area  (in^2) 
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%%%% 

Lmaxs 
WmaxB 
CgL  = 


Constraint 

Parameters 

72; 

%  max  total  length 

30; 

%  max  total  weight 

=  1/3; 

%pod  forward  Cg  limit  (frac  of 

lug  spacing  back  from  front  lug) 


%%%%  Rocket  Parameters 

Lr  »  55.125; 
CgR_live  «  29.96; 
CgR__burn  =  33.55; 
WR_live  =  22.95; 

WR__burn  =  15.73; 


%  rocket  length 

%  live  rocket  Cg  (inches  from  base) 
%  fired  rocket  Cg  (inches  from  base) 

%  live  rocket  mass  (Ibm) 

%  fired  rocket  mass  (Ibm) 


%%%%  Launcher  Pod  Parameters 

Llug7  s  43.2; 

Llugl9=  43.4; 

Cg7  =  Llug7-(CgL)*14; 
Cgl9sLlugl9- (CgL) *14; 
Cgpod  e  32.6; 

Wpod7  =  196.2; 

Wpodl9=  516.0; 


%  front  lug  from  back  of  launcher-?  (in) 

%  front  lug  from  back  of  launcher— 19  (in) 

%  constraint  cg  (e.g.  within  middle  half  of  lug  space) 
%  constraint  cg  (e.g.  within  middle  half  of  lug  space) 
%  loaded  pod  Cg  from  back  w/o  kits  (in) 

%  loaded  pod (7)  mass  w/o  kits  (Ibm) 

%  loaded  pod (19)  mass  w/o  kits  (Ibm) 


%  (VARIOUS  DENSITY  CURVES  PLOTTED) 


Wk  =  Wkl  *  *ones  ( 1 ,  length  (Dkl ) )  ; 
Lk  =  Wk(:,l)*(l./(area*Dkl))  ; 
Cgk=  Lr+(l-Cgkl)*Lk; 


%  matrix  of  kit  masses 
%  matrix  of  kit  lengths  (in) 

%  kit  Cg  locations  (inches  from  base) 


Mk  =  Wk.*Cgk; 


%  Moment  about  base  due  to  kit  (Ibm  in) 


%  centers  of  gravity  are  calculated  from  base  in  inches 
Cg_live  =  (WR_live*CgR_live+Mk)  .  /  (WR_live+Wk)  ; 

Cg_bum  =  (WR_bum*CgR_bum+Mk)  .  /  (WRjDum+Wk)  ; 


%  centers  of  gravity  are  changed  to  calibers  from  nose 
Cg_l ive  =  ( ( Lr +Lk )  -Cg__l ive )  / dia  ; 

Cg_bum  =  ( (Lr+Lk) -Cg_bum) /dia; 

%figure 
figure  (1), 

plot  (Wk( : ,  1) ,  Cg_live)  ,  grid  on 
title  ( '  Live  Round  with  Kit ' ) 
ylabel  ('Cg  (Calibers  from  nose)  ' ) 
xlabel  ('kit  weight  (Ibm)') 
hold  on 


69 


figure  (2), 

plot  {Wk(:,l),  Cg_bum) ,  grid  on 
title  ( '  Fired  Round  with  Kit ' ) 
ylabel  {'Cg  (Calibers  from  nose)  ' ) 
xlabel  ('kit  weight  (Ibm)  ' ) 
hold  on 


%  (VARIOUS  LENGTH  CURVES  PLOTTED) 
Lk  =  Lkl; 

Dmin  =  0.8*min(Dkl)  ; 

Dmax  =  1.2*max(Dkl)  ; 

Dk  =  linspace (Dmin,  Dmax)  ; 

Dk  =  Dk ' *ones ( 1 , length (Lk) ) ; 

Wk  =  area*Dk( : , 1) *Lk; 

Cgk=  Lr+(l-Cgkl)*Lk; 

Cgk=  ones ( length ( Wk ) , 1 ) *Cgk ; 

Mk  =  Wk,*Cgk; 


%  specific  kit  lengths  (in) 

%  lower  limit  on  assumed  density  range 
%  upper  limit  on  assumed  density  range 
%  assummed  kit  density  range  (ltm/in^3) 
%  matrix  of  kit  densities 
%  matrix  of  kit  masses  (Ibm) 

%  kit  Cg  locations  (inches  from  base) 

%  matrix  of  Cg  locations 

%  Moment  about  base  due  to  kit  (Ibrn  in) 


%  centers  of  gravity  are  calculated  from  base  in  inches 
Cg_live  =  (WR_live*CgR_live+Mk)  .  /  {WR_live+Wk)  ; 

Cg_jDum  =  (WRjDum*CgR_bum+Mk)  .  /  (WR_bum+Wk)  ; 


%  centers  of  gravity  are  changed  to  calibers  from  nose 
Cg_live  =  ( (Lr+ones  (length (Wk) ,  1)  *Lk) -Cg_live) /dia; 

Cg_bum  =  ( (Lr+ones  (length  (Wk) ,  1)  *Lk)  -Cg_bum)  /dia; 

figure  (1), 

plot (Wk, Cg_live) 

plot  ([  (Wnnax-WR_live)  *ones(  1,2)  ] ,  [1,20])  %itBx  weight  constraint 
figure  (2), 

plot  (Wk,  Cg_Jxim) 

plot  ( [  (V*nax-WR_live)  *ones  (1,2)],  [1,20])  %max  weight  constraint 
%  (POD  CENTER  OF  GRAVITY  CONSTRAINT) 

Wk  =  Wkl '  *ones  (1, 2) ;  %  matrix  of  kit  masses 

%length  from  loaded  pod  (w/o  kits)  cg  to  constraint  Cg  (in) 

LI  =  ones(length(Wk),l)*([Cg7,Cgl9]“Cgpod); 

%  loaded  pod  weights  (Ibm) 

Wp  =  ones ( length ( Wk)  ,l)*[Wt>od7,V^odl9]  ; 

X  =  ones  (length (Wk) ,  1)  *  [7 , 19] ;  %  number  of  kits  per  pod 

%  length  from  rocket  nose  back  to  front  lug 
%  (assume  rocket  base  4"  from  back  of  pod) 

Lt  =  ones (length (Wk) , 1) * [Lr-Llug7+4,Lr-Llugl9+4] ; 
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%  Solve  for  allowable  kit  length 

Lk  =  (l/(l-Cgkl))*((Ll.*Wp)./(X.*Wk)-Lt); 

Cgk=  Lr+ (l“Cgkl) *Lk;  %  kit  Cg  locations  (inches  from  base) 

%  center  of  gravity  calculated  from  base  in  inches 
Cg_live  =  {  (WR_live*CgR_live)  +  (Cgk.  *Wk)  )  .  /  (WR_live+Wk)  ; 

Cg_burn  =  (  (WR_burn*CgR_bum)  +  (Cgk.  *Wk) )  .  /  (WR_burn+Wk)  ; 

%  center  of  gravity  changed  to  calibers  from  nose 
Cg_live  =  ( (Lr+Lk) -Cg^live) /dia; 

Cg_burn  =  ( (Lr+Lk) ~Cg_burn) /dia; 

figure  (1) , 

plot (Wk, Cg_live) 
axis{11.8,7.2,7.8,14]) 

%hold  off 
figure  (2)  , 

plot (Wk, Cg_burn) 
axis([1.8,7.2,6.8,13]) 

%hold  off 

E.  STEADY-TURN  MANEUVER 

%  This  m-file  computes  thrust  required,  T,  for  certain  load  factor,  n 
%  (g's),  in  a  steady  turn  for  various  Mach  niimbers.  The  load  factor  is 
%  achieved  using  a  thruster  forward  of  the  CG. 

% 

%  The  inputs  include: 

%  Scenario  Definition: 

%  n:  desired  load  factor  (g's) 

%  rt:  desired  turn  radius  (ft) 

%  alt:  altitude  (ft)  (std  day  assumed) 

%  Mnin:  min  Mach  number 

%  Mnax:  max  Mach  number 

%  Rocket  Configuration: 

%  Wght:  rocket  mass  (Ibm) 

%  Dia:  rocket  cross-section  diameter  (in)  (for  ref  area  and  length) 

%  Xcg:  rocket  center  of  gravity  in  calibers  from  nose 

%  Xt:  thruster  location  in  calibers  from  nose 

%  Aerodynamic  Parameters : 

%  CNa:  normal  force  coefficient  per  radian  of  alpha  (f(Mach)) 

%  Xcp:  center  of  pressure  location  in  calibers  from  nose  [f (Mach) 3 

%  CMq:  moment  coefficient  per  rad/sec  of  pitch  rate  [f(Mach)] 

%  (Aero  Paramaters  are  linearly  interpolated  between  points  in  table) 
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%  Code  uses  normal  force  and  pitch  moment  equations: 

%  Fz=nW  =  N  +  T=  (Qbar)  *S*CNa*alpha  +  T,  and 
%  M  =  0  =  -N*(Xcp-Xcg)  +  T*  (Xcg-Xt)  -  CMq*q*  (Dia/2V)  *Qbar*S*Dia 
%  {since,  q  =  theta-dot  =  V/R  =  V/(V^2/n)  =  n/V 

%  and  N  =  (Qbar)  *S*CNa*alpha} 

%  M  =(Qbar*S*CNa*alpha)  (Xcp-Xcg)  +T*  (Xcg-Xt)  -  CMq*n*(Dia^2/V^2)*Qbar*S 
%  The  two  equations  are  solved  for  the  unknowns,  alpha  and  T,  using  several 
%  values  of  Mach  number. 


%  W.  Mark  Wonnacott,  Naval  Postgraduate  School,  May  1997 
clear,  close  all 


%  INPUTS 
%  ====== 

%  Scenario  Definition 

n  =  15; 

Rt  ~  6000; 
req_type  =  2 ; 
alt  B  500; 

MmixiB  0.4; 

HmaxB  1.6; 


%  desired  load  factor  (g's) 

%  desired  min  turn  radius  (ft) 

%  requirement  (l:tum  rate,  2: turn  radius) 
%  std  day  pressure  altitude  (ft  MSL) 

%  minimum  Mach  number  for  calculations 
%  maximum  Mach  number  for  calculations 
%  Rocket  Configuration  (based  on  bumed-out,  30  Itam,  68.75"  rocket) 

Wght  =  22.78;  %  rocket  mass  (Ibm) 

Dia  =  2.75;  %  reference  length  (in) 

Xcg  =  10.0;  %CGin  calibers  from  nose 

Xt  =  2.2;  %  thruster  location  in  calibers  from  nose 

%  Aerodynamic  Parameters  (based  on  final  config  68.75"-w/  3  case#8  strakes) 
% 

aero_tabla  = 


Mach# 

CNa 

Xcp 

CMq 

to 

15.0616 

11.52547 

-4166.814 

0.6614 

16.41589 

11.32152 

-4556.4946 

0.866 

18.10928 

11.01214 

-5041.0303 

0.9682 

20.04303 

10.36514 

-5518.2495 

1 

20.05156 

10.35951 

-5519.5352 

1.0308 

20.05156 

10.35951 

-5519.5352 

1.118 

18.91632 

10.45431 

-4859.3901 

1.25 

17.51701 

10.5577 

-4376.4146 

1.4142 

16.10659 

10.68451 

-3938.8662 

1.6008 

14.68889 

10.88046 

-3541.0081 

1.8028 

13.47025 

10.99184 

-3202.1553 

2 . 0156 

12.38573 

11.05396 

-2916.2678] 

Mach_dat=  aero_table  ( : ,  1 )  ; 
CNa_dat  =  aero_table  ( : ,  2 )  ; 
Xcp_jiat  =  1.03  *aero_table  ( : ,  3 ) ; 
CMq_dat  =  aero_table  ( : ,  4 )  ; 


%  Mach  numbers 

%  Bo<^  Normal  Force  per  Radian  Alpha 
%  Center  of  Pressure  in  calibers  from  nose 
%  Pitching  Moment  per  rad/sec  of  Pitch  Rate 
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%  CONSTANTS 
gO  =  32.17; 

TenpO  =  518.67; 
PresO  =  14.696*144; 
B_atm=  0.003566; 
R_air  =  53.35; 
k_air  =  1.4; 


%  Gravity  Constant. 

%surface  tenperature  (°R) 
%surface  pressure  (lbf/ft^2) 
%teirp  lapse  rate  (®R/ft) 

%gas  constant  (ft  lb/(lbm  °R) ) 
% specific  heat  ratio 


%  INTERMEDIATE  CALCULATIONS 

%  Speed  of  Sound 

tenp  =  TeitpO-B_atm*alt; 

pres  =  PresO*  (1-B_atm*alt/Tenp0)  ^  (1/R_air/B_atin)  ; 
rho  =  pres/ (R_air*tertp)  ; 


Vs  =  sqrt  (k_air*R_air*terrp*gO) ; 

%  Rocket  Config 

d  =  Dia/12;  % 
S  =  pi*d'"2/4;  % 
XT  =  (Xcg  -Xt)*d;  % 
XP_dat  =  (Xcp_dat“Xcg)  *d;  % 
Fz  =  n  *  Wght;  % 
%  Mach  nuinbers  and  speeds 

%  adjust  min  and  max  Mach  to  table  values 
Mmin  =  max{min(Mach_dat)  ,Mnain) ;  % 
Itoiax  =  min{max{Mach_dat)  ,t&nax) ;  % 
kmax  =50;  % 
Mach  =  1  inspace  (Mnin^Jfinax,  kmax) ;  % 


%  terrperature  at  altitude  ('^) 
%pressure  @  alt  (lb/ft^2) 

%  density  at  altitude  (lbm/ft^3) 

%  speed  of  sound  (ft/sec) 

reference  length  (ft) 

reference  area  (ft^2) 

thruster  location  from  CG  (ft) 

center  of  pressure  location  from  CG  (ft) 

normal  force  (lb) 

min  Mach  number  to  be  used 
max  Mach  number  to  be  used 
number  of  Mach  nuinbers  to  be  used 
vector  of  Mach  numbers  to  be  used 


%  THRUST  AND  ALPHA  CALCULATIONS  FOR  EACH  MACH  NUMBER 


for  k  =  l:kmax 
V  =  Mach(k)*Vs; 

QS  =  S  *  l/(2*gO)*rho*V"2; 

CNa  =  interpl(Mach_dat,CNa_dat,Mach(k) )  ; 
Xcp  =  interpl(Mach_dat,XP_dat,Mach(k) )  ; 
CMq  =  interpl(Mach_dat,CM^dat,Mach(k) )  ; 
A  =  [QS*CNa,  1; 

-QS*CNa*Xcp,  XT]; 
if  req_type  ==  2, 

B=  [(V"2*Wght)/(gO*Rt); 

(1/2) *QS*CMq*d-2/Rt]  ; 

else 

B=  [Fz;  (l/2)*QS*CMq*n*gO*(d/V)^2]  ; 
end 

X  =  inv(A)*B; 
a(k)  =  X(l)*180/pi; 

T(k)  =  X(2); 
end 


%  rocket  speed  (ft/s) 

%  (dyn  pres— lb/ ft ^2)  (ref  area— ft ^2)  (lb) 
%CNa  for  Mach# 

%Xcp  for  Mach# 

%CMq  for  Mach# 


%  RHS  vector  for  turn  radius 


%  RHS  vector  for  turn  rate 

%  solution  vector 
%  alpha  in  degrees 
%  required  thrust  in  pounds 
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subplot  (2,1,1) 

plot  (Mach,T, 'r ' ) ,  grid  on 
xlabel  ( '  Mach  Number ' ) 
ylabel(  'Required  Thrust  (Ibf)  ’ ) 
hold  on 

subplot  (2,1,2) 

plot  (Mach, a, 'r ' ) ,  grid  on 
xlabel  ( '  Mach  Number ' ) 
ylabel (' Required  alpha  (deg)’) 
hold  on 


F.  SEEKER  FIELD  OF  VIEW 

%  This  m^file  calculates  the  probability  of  a  target  being  in  the 
%  seeker  field  of  view,  FOV  (deg) ,  for  for  several  seeker  acquis t ion 
%  ranges  (Racq)  and  fields  of  view.  Only  single  launch  range  is  used. 

% 

%  inputs  include:  Launcher  Pointing  Error,  Pnt  (deg)  -  1  sigma 
%  Missile  unguided  ballistic  error,  B  (mils)  -  1  sigma 

%  Missile  distance  time  profile,  (km)  (function  of  fly  out  time) 

%  Guidance  time  delay,  Gd  (sec) 

%  Time  delay  for  target  ID,  td  (sec) 

%  Min  and  Max  Target  ^eed  (mph) 

%  also  required  are: 

%  Seeker  acquis tion  ranges  (km)  -  (vector) 

%  Min  and  Max  FOV  (deg) 

%  Launch  Range  ( km) 

%  Mark  Wonnacott,  April  1997 
%  Naval  Postgraduate  School 
clear,  close  all 

%%  Target  Description 

VTmin  =  20;  %  (nph)  min  target  speed 

VTmax  =  60;  %  (mph)  max  target  speed 

%%  Missile  Configuration  Inputs 
%  Launcher  Pointing  Error 

Piib  =  3;  %  (deg~  1  sig)  initial  pointing  error 

%  Unguided  Missile  Ballistic  Error 

B  -  12;  %  (mills-  1  sig)  meters  error  per  km  flight 

%  Guidance  Delay 

=1.0;  %  (sec)  max  range  delay  (no  guidance  for  1  sec) 

%  Delay  time  to  ID  target 

td  =  5*1/8;  %  (sec)  5  laser  pulses  received  at  8  Hz 
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%  Seeker  Field  of  View 
FOVmin  =  .5; 

FOVuax  =  45; 

%  Seeker  Acquisition  Range 
Racq  =  [1.15,1.5,2,3,6] ; 

%  Launch  Range 
RL  =  6; 

%  Missile  Speed  Table 


%  (deg)  min  half  FoV 
%  (deg)  max  half  FoV 

%  (km)  parametric  seeker  acquis t ion  ranges 
%  (km)  maximum  launch  range 


Time  -  [0,0.22,0.89,1.02,1.67,...  %f light  time  (s) 

2.52,4.96,7.91,11.24,14.88,18.8,22.99,27.3] ; 

Dist  =  [0,0.01,0.2,0.27,0.6,1,2,3,4,5,6,7,8];  % flight  dist  (m) 


%%  Inputs  converted  to  usable  variables 
VTmin2  =  VTmin/2237; 

VTmax2  =  VTmax/2237; 

PntR  =  Pnt*pi/180; 

Bkm  =  B/1000; 

kmax  =  length  (Racq)  ; 

raraax  =  30; 

FOV  =  linspace(FO\^n,FO\^nax,mmax)  *pi/180; 
Rg  =  interpl  (Time,  Dist,  Gd); 


%  min  target  speed  (mph->km/s) 

%  max  target  speed  (mph->km/s) 

%  initial  pointing  error  (deg“>rad) 
%  X-trk  error  (km  per  km  flight) 

%  number  of  seeker  acq.  ranges 
%  number  of  FOV's 
%  FOV ’ s  ( deg- >rad ) 

%  unguided  flyout  range  (km) 


Pacq  =  zeros  (kmax, mmax) ; 


%  pre-set  Pacq  matrix  size 


for  k  =  l:kmax 
for  m  =  limmax 

D  =  Racq(k)  *cos  (FOV(m) )  ; 
F  =  Racq  (k)*sin  (FOV (m) )  ; 
if  RL-Rg  >  D 
Rd  =  RL  -  D; 
else 

Rd  =  RL  -  Rg; 
end 


%  calculations  for  each  Racq 
%  calculations  for  each  FOV 

%  detection  range  (trgt  ©  FOV  limit) 

%  max  x-trk  offset  for  FOV 

%  range  flown  before  detection 

%  range  flown  when  guided  flight  begins 


tD  =  interpl  (Dist, Time,  Rd)  ; 
tT  =  tD  +  td; 

Rt  =  interpl  (Time, Dist,  tT) ; 
s  =  Rt*sqrt(tan(PntR)^2+(Bkm)^2) 
rl  =  tT  *  VTmin2; 
r2  =  tT  *  VTraax2; 


%  time  at  detection 
%  time  at  start  of  turn 
%  distance  flown  at  start  of  turn 
%1  sig  X-trk  error  at  start  of  turn 
%  min  target  motion  at  start  of  turn 
%  max  target  motion  at  start  of  turn 


L5  =  min(5*s,  r2+F) ;  %  integration  upper  limit 

%%%  Integration  of  Pr  (Xr)  *Pt  (Xr-F,Xr+F) 

%%%  from  0  to  L5  done  using  'quad8'  (Pacq  function  in  'ProhAcq') 
Pacq(k,m)  =  2*quad8  ( 'ProhAcq' ,0,L5, 0.01,  []  ,F,s,rl,r2)  ; 
end 
end 
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FOV  =  FOV*180/pi;  %  field  of  view  changed  back  to  degrees 

figure 

plot  (FOV, Pacq) ,  grid  on 

title  (['Vt  =  [ '  ,nuin2str  (VTmin) ,  *  ->  '  ,num2str  (VTmax) ,  ’  ]  itph']) 

legend (nuin2str{Racq(l) )  ,nuin2str  (Racq(2) )  ,num2str  (Racq(3) )  ,nuin2str  (Racq(4) )  ,num2str  (Rac 
q(5))) 

ylabel  ('Probability  of  Teirget  in  FOV), 
xlabelChalf  FOV  (±deg)  ' ) 


ProbAcq  Function 

fionction  Pacq  =  ProbAcq (Xr, F, s,rl,r2) 


%  Returns  the  product  of  the  assumed  probability  of  the  rocket  being 
%  approximately  at  a  certain  cross- track  position,  Xr,  and  the 
%  probability  of  the  target  being  between  ±F  to  either  side  of 
%  the  rocket  position  (seeker's  FOV  limits)  . 

%  The  probabilities  are  based  on  the  1-sigma  rocket  cross-track  error— s 
%  and  the  min  and  max  possible  target  cross-track  motion— rl  and  r2. 

% 

%  This  function  has  to  take  Xr  as  a  vector  and  output  Pacq  as  a  vector!  !  ! 

%  W.  Mark  Wonnacott,  Naval  Postgraduate  School,  1997 

if  rl  ==  0;  rl  =  0.0001;  end  %  avoid  division  by  zero 

Prkt  =  1/ (s*sqrt  (2*pi) )  *exp(-Xr.^2.  /  (2*s^2) ) ;  %normal  distribution  of  rocket 


%  Find  probability  of  target  in  FOV  extremes  (Pt(xl,x2))  for  each  Xr 
for  k  =  1: length (Xr) 


X  =  Xr(k) ; 

XI  =  X  -  F; 

X2  =  X  +  F; 
if  XI  >  r2 

Ptgt(k)  =  0; 
elseif  XI  <  -r2 
Ptgt(k)  =  1; 
else 


%  step  through  each  rocket  x-trk  position 
%  extreme  bound  of  FOV 
%  extreme  bound  of  FOV 

%  weapon  farther  out  than  possible  for  target 
%  all  possible  target  locations  inside  FOV 


%  conpute  area  bounded  rl  &  r2  circles,  fly-out  axis,  and  X2 
LI  =  min(X2,r2);  %  limit  of  intgrt'n  for  r2 

L2  =  min(X2,rl);  %  limit  of  intgrt’n  for  rl 

Ala  =  Ll*sqrt  (r2^2-Ll'^2)  +  r2^2*asin(Ll/r2) ;  %bounded  area  in  r2  circle 

Alb  =  L2*sqrt  (rl^2-L2^2)  +  rl^2*asin(L2/rl) ;  %bounded  area  in  rl  circle 


circle 

circle 


A1  =  Ala-Alb; 
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%  conpute  area  bounded  rl  &  r2  circles,  fly-out  axis,  and  XI 
L3  =  niin(abs  (XI)  ,r2)  ; 

L4  =  inin(abs(Xl)  ,rl)  ; 

A2a  =  L3*sqrt{r2'^2-L3^2)  +  r2'^2*asin(L3/r2)  ; 

A2b  =  L4*sqrt(rl"2-L4"^2)  +  rl^2*asin(L4/rl)  ; 

A2  =  A2a-A2b; 

A  =  Al-sign(Xl)  *A2; 

Ptgt{k)  =  A  /  (pi*(r2^2-rl"2)) ; 

end 

end 

Pacq  =  Prkt.*Ptgt*100; 
return,  end 

G.  6DOF  CODES 

1.  Simulation  Initilization 

%  File  Initializes  Simulation 
clear 

data_out  =  [2001,5]; 

%%%  Missile  INITIAL  CONDITIONS  (6  DOF  Flight) 

%%%  (u ,  V,  w,  p ,  q,  r ,  phi ,  theta ,  psi ,  p_rocket ,  pn,  pe ,  ph) 

%  ground  launch  initial  conditions 

Msl_Xo  =  [0;0;0;0;0;0;0;10*pi/180;0;0;0;0;0]  ; 

%  target  pointing  after  12  sec  flight  (10°  elev  launch)  trgt  @  5000m 
%Msl_Xo  =  [285;0;0;0.03; .01;0;0;-.383;0;-125.6;4230;0;310] ;%theta  =0.0246 
%  straight  and  level  cruise  after  12  sec  flight 
%Msl_Xo  =  [300;0;0;0;0;0;0;0;0;-125.6;0;0;300]  ; 

%%%  Target  INITIAL  CONDITIONS  (6  DOF  Flight) 

%%%  (Pn,  Pe) 

Tgt_Xo  =  [5000; 001; 

T_Vn  =  0; 

T_Ve  =  0; 

Talt  =  0; 


Aerodat  %  weight  and  aero  data  file  for  2.75"  HRGK -equipped  rocket 
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2.  Aerodat 

%  This  File  Contains  Look-up  Table  Data 
%  2.75"  Rocket  with  Final  Design  Guidance  Kit 

global  gO  TerrpO  PresO  B_atin  R_air  S  d  Sd  Xt  mu  rho  c; 


%CONSTANTS  &  CONVERSIONS 
%%%  Conversions 
in2m  =  39.37; 
lb2kg  =  2.2046; 

%%%  Gravity 
gO  =  9.81; 

%%%  Atmospheric  Data: 
TerrpO  =  288.15; 

PresO  =  101330; 

B_atm=  0.00650; 

R_air  =  287.0; 


%  inches  per  meter 
%pounds  per  kilogram 

%  Gravity  Constant,  g'O,  in  m/s^2. 

%surface  tenperature  (K) 

%surface  pressure  (Pa) 

%terrp  lapse  rate  {°C/m) 

%gas  constant  (N  m/(kg  K) ) 


%  Guidance  Package  Pcirameters 

radius  =  2.75/2/in2m;  %Guidance  Package  Radius  (m) 


GP_Mass0  =  7/lb2kg; 
GP_Massf  =  7/lb2kg; 
GP_Lngth  =  13.625/in2m; 


%  Initial  Guidance  Package  Mass  (kg) 
%Final  Guidance  Package  Mass  (kg) 
%Guidance  Package  Length  (m) 


Xt  =  6.6/in2m; 


%Thruster  Location  from  Nose  (m) 


mu  =  0.006; 


%Rocket  to  Kit  Coefficient  of  Friction 


GP_CG0  =  8.5/in2m;  %Initial  Guidance  Package  CG  from  nose  (m) 

GP_CGf  =  8.5/in2m;  %Final  Guidance  Package  CG  from  nose  (m) 


GP_Jxx0  =  ...  %Aprox. Initial  Axial  Moment  of  Inertia  (kg  m^2) 

GP__Mass0*radius^2  /2  ; 

GP_Jxxf  =  —  %;^rox. Final  Axial  Moment  of  Inertia  (kg  m^2) 

GPJMassf  *radius'"2  /2  ; 


GP_jJyy0  =  —  %^rox. Initial  Transverse  M.  of  I.  (kg  m^2) 

GP__Mass0*(GP_CG0--GP_Lngth/2)'^2  +  GP_MassO*GP_Lngth^2/12; 
GP_jJyyf  =  ...  %i^rox. Final  Transverse  M.  of  I.  (kg  m^2) 

GPjyiassf*(GP_CGf-GP_Lngth/2)^2  +  GP_Massf  *GP_Lngth^2/12  ; 
FlyTime  =  30;  %i^proximate  Flyout  Time  Bound  (sec) 
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%  Rocket  Parameters  from  IHSP  89-289,  Mk-66  Data  Book 

%Rocket  Diameter  (m) 


d  =  2.75/in2m; 

R_MassO  =  22.95/lb2kg; 

RJiassf  =  15.73/lb2kg; 

R^Lngth  =  55.125/in2m; 

R_CG0  =  ... 

R— Lngth-29 . 96/in2m  +  GP_Ijngth; 
R__CX5f  =  ... 

R_Lngth-33 .55/in2m  +  GP_Lngth; 


%Initial  Rocket  Mass  (kg) 

%Final  Rocket  Mass  (kg) 

%Rocket  Length  (m) 

% Initial  Rocket  CG  from  nose  (m) 

%Final  Rocket  CG  from  nose  (m) 


R_JxxO  =  26.2/(lb2kg*in2m^2); 
R_Jxxf  =  19.7/(lb2kg*in2m^2)  ; 


%Initial  Axial  Moment  of  Inertia  (kg  m'^2) 
%Final  Axial  Moment  of  Inertia  (kg  m^2) 


R_JyyO  =  6248/(lb2kg*in2m^2) ;  %Initial  Transverse  M.  of  I.  (kg  m^2) 

R— Jyyf  =  5008/ (Ib2kg*in2m^2) ;  %Final  Transverse  M.  of  I.  (kg  m^2) 


BumTime  =  1.05; 

%Average  Rocket  Motor 

Bum  Time 

(sec) 

ThrustHot_j3at= . . . 

%Thrust  Time  History 

(sec  &  Ibf 

for  150'^F) 

[0.00  0 

0.40 

1520 

0.85 

1800 

0.05 

1600 

0.45 

1560 

0.90 

1700 

0.07 

1700 

0.50 

1610 

0.95 

600 

0.10 

1690 

0.55 

1650 

1.00 

200 

0.15 

1600 

0.60 

1700 

1.05 

40 

0.20 

1530 

0.65 

1730 

1.10 

10 

0.25 

1500 

0.70 

1750 

1.15 

0 

0.30 

1490 

0.75 

1780 

FlyTime  0]; 

0.35 

1510 

0.80 

1800 

ThrustCold_dat=. . . 

%Thrust  Time  History  (sec  &  Ibf  for  -50°F) 

[0.00  0 

0.4 

1380 

0.85 

1680 

0.05 

1150 

0.45 

1350 

0.9 

1760 

0.07 

1300 

0.5 

1340 

0.95- 

1750 

0.1 

1380 

0.55 

1350 

1.00 

1600 

0.15 

1400 

0.6 

1350 

1.05 

700 

0.2 

1400 

0.65 

1400 

1.1 

100 

0.25 

1400 

0.7 

1440 

1.15 

0 

0.3 

1390 

0.75 

1500 

FlyTime  0] ; 

0.35 

1400 

0.8 

1590 

ThrustHot_dat  ( : ,  2 )  =  ThrustHot^dat  ( : ,  2 )  *g0/lb2kg; 
ThrustCold_dat( :  ,2)  =  ThrustCold_dat  ( : , 2)  *g0/lb2kg; 
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%  Reference  Area/Length  and  Mass  &  Moment  Time  Histories 

S  =  pi*(d^2)/4;  %  Cross-Sectional  Area  (Ref  Area-^'^2) 

Sd  =  S*d;  %  Ref  Area  x  Ref  Length  {m^3) 

t  =  [0;  BumTime;  FlyTime] ;  %time  vector  for  time  histories 

Mass_dat  =  [t,  [RJMassO+GP_MassO  %  Mass  Time  History  (sec  &  kg) 

R__Mas  s  f +GP_Mas  s  0 
R_Massf+GP_Massf  ]  ]  ; 


cg0=  (R_CG0*R_Mass0+GP_CG0*GP_Mass0)  /  (RJMassO+GP_MassO)  ; 
cgl=  (R_CGf  *RJMassf+GP_CGO*GP_MassO)  /  (R_Massf+GP_MassO)  ; 
cg2=  (R_CGf *R_Massf+GP_CGf *GP_Massf )  /  (R_Massf+GP_Massf )  ; 


CG_dat  =  [t, 


[cgO 

cgl 

cg2]]; 


%  CG  History  (sec  &  m  from  nose) 


JxO  =  R_JxxO+GP_JxxO; 

Jxl  =  R_jJxxf +GP_JxxO  ; 

Jx2  =  R__Jxxf +GP_Jxxf  ; 

Jx„dat  =  [t,  [JxO  %  Axial  M.  of  I.  History  (kg  m^2) 

Jxl 
Jx2]]; 

JxR_dat  =  [t,  [R__JxxO;  R_Jxxf;  R_Jxxf  ]  ]  ; 

JxK_dat  =  [t,  [GP_JxxO;GP_JxxO;GP_Jxxf]]; 


JyO  =  R_Jyy0+R_Mass0*(R_CG0-cg0)^2+GP_Jyy0+GP__Mass0*(GP_CG0-cg0)^2; 
Jyl  =  R_Jyyf+R_Massf  *  (R_CGf-cgl)  ^2+GP_JyyO+GP_MassO*  (GP_CG0-cgl)  ^2  ; 
Jy2  =  R_Jyyf +R_Massf  *  (R_CGf-cg2 )  ^2+GP_Jyyf +GP__Massf  *  (GP_CGf -cg2 )  ^2  ; 
Jy_dat  =  [t,  [JyO  %  M.  of  I.  History  (kg  m^2) 

Jyl 

Jy2]]; 

%  Rocket  Roll  Acceleration  History 


R_Pdot 

=  . . . 

% 

Rocket  Roll  Accel  (rev/s/s) 

[0 

0 

1.8 

-8.0 

5.8 

-4.5 

0.5 

45.0 

2.0 

0.0 

6.2 

-12.5 

0.6 

100.0 

2.4 

2.0 

6.6 

1 

O 

O 

0.9 

15.0 

3.0 

7.5 

7.0 

-7.0 

1.0 

-50.0 

3.6 

8.5 

8.0 

-4.0 

1.2 

-120.0 

4.0 

9.5 

9.5 

0.0 

1.4 

-92.0 

4.4 

7.5 

12. 

1.0 

1.6 

-34.0 

5.6 

0.0 

21. 

0.1] 

R_Pdot ( 

:,2)  =  R_Pdot(:,2)*2*pi; 

% 

converted  to  rad/s/s 

80 


%  Aerodynamics  Data  from  ADAM'S  MADR 


Mach_dat  =  ... 


[0  1 

0.6614  1 

0.8660  1 

0.9682  1 


.0000 

1.4142 

.0308 

1.6008 

.1180 

1.8028 

.2500 

2.0156] 

%Static  Coefficient/Derivative  Components  Based  on  Mach  Number  Vector 


Cd0Coast_dat  = 
[0.52292 
0.48521 
0.51417 
0.75708 

Cd0Bum_dat  =  ... 
[0.42291 
0.38521 
0.41417 
0.65708 
CNa__dat  =  ... 
[15.0616 
16.41589 
18.10928 
20.04303 
CP_dat  =  . . . 
[11.52547 
11.32152 
11.01214 
10.36514 

CP_dat  =  CP_dat*d; 


%  Skin,  Wave,  and  Base  Drag 
0.83218 
1.0067 
1.00424 
1.02549 


(rocket  coasting) 

1.09629 

1.13922 

1.16393 

1.17799]; 


%  Skin  and  Wave  Drag  (rocket  burning) 
0.73219 
0.81307 
0.81061 
0.83187 

%  Body  Normal  Force  per  Radian  Alpha  or  Beta 
20.05156 
20.05156 
18.91632 
17.51701 

%  Center  of  Pressure  in  Calibers 


0.90267 

0.95226 

0.99252 

1.02295]; 

16.10659 

14.68889 

13.47025 

12.38573]; 


10.35951 

10.35951 

10.45431 

10.5577 


10.68451 

10.88046 

10.99184 

11.05396]; 


%pyanamic  Coefficients /Derivatives  Based  on  Mach  Number  Vector 
% - 


CMqudat  =  ... 

[-4166.814 
-4556.4946 
-5041.0303 
-5518.2495 
Clp_dat  =  , .  . 

[-7.513 

-7.513 

-7.513 

-7.513 


%  Pitching  damping  per  Radian/ sec  of  Pitch  Rate 


-5519.5352 

-5519.5352 

-4859.3901 

-4376.4146 


-3938.8662 

-3541.0081 

-3202.1553 

-2916.2678] 


%  Roll  damping  (HRGK  only)  /Radian/sec  of  Roll  Rate 


-7.513 

-7.513 

-7.2767 

-7.0348 


-6.7529 

-6.453 

-5.5802 

-4.7831]; 
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3.  Equations  of  Motion 

function  xo  =  MotionEqns  (uo) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  This  Function  coirputes  the  Flat-Earth,  Body  Axes, 

%%%  6  “DOF  Dynamics  Equations  for  a 

%%%  Thruster “Controlled  Missile 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  Adapted  for  rolling  airframe,  thrust  control,  and 

%%%  variable  aero  coefficient  and  missile  parameters  by 

%%%  W.  Mark  Wonnacott,  M^ch  1997 

%%%  Naval  Postgraduate  School 

%%%  Adapted  from  the  G.  Hutchins'  Code  for  EC  4330. 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  Original  Code  hy  R.  G.  Hutchins,  10  December  1996 
%%%  Thanks  to 

%%%  Rob  King  and  Mark  Wonnacott  for  aero  assistance 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  The  State  Vector  is  defined  as: 

%%%  Body“Axes  Velocities 


%%% 

x(l)  =  U, 

speed  along  missile  axis,  x  direction 

%%% 

x(2)  =  V, 

speed  of  sideslip,  y  direction 

%%% 

x(3)  =  W, 

speed  along  z  direction  ("Down") 

%%% 

Body“Axes  Angular  Velocities 

%%% 

x(4)  =  P, 

roll  rate 

%%% 

x(5)  =  Q, 

pitch  rate 

%%% 

x(6)  =  R, 

yaw  rate 

%%% 

Euler  Angles 

%%% 

x(7)  =  phi, 

,  ’  roll  relative  to  vertical  down 

%%%  x(8)  =  theta,  pitch 

%%%  x(9)  =  psi,  yaw 

%%%  Rocket  Rolling 

%%%  x(10)  =  Pr,  rocket  roll  rate 

%%%  Location  Variables 

%%%  x(ll)  =  Pn,  Position  North  of  (0,0,0) 

%%%  x(12)  =  Pe,  Position  East  of  (0,0,0) 

%%%  x(13)  =  h.  Height 


%%%  The 

input 

vector 

is  defined  as: 

%%% 

U(l) 

=  Tl, 

thrust  from  thruster  #1 

%%% 

u(2) 

=  T2, 

thrust  from  thruster  #2 

%%% 

u(3) 

=  T3, 

thrust  from  thruster  #3 

%%% 

u(4) 

=  T, 

thrust  through  the  missile  x  axis 
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%%%Other  Missile  Parameter  Inputs  (fimctions  of  Mach  or  time)  : 
%%%  P{1)  =  mass,  missile  mass  (kg)  [f(t)] 

%%%  P{2)  =  Xcg,  missile  Cg  (m  from  nose)  [f(t)] 

%%%  P(3)  =  Xcp,  center  of  pressure  (m  from  nose)  [f{M)] 

%%%  P(4)  =  JxR  rocket  MofI  about  long  axis  (kg  m'^2) 

%%%  P(5)  =  JxK  kit  MofI  about  long  axis  (kg  m'^2) 

%%%  P(6)  =  Jy=:Jz,  missile  Mofl-y  or  z  axes  [f(t)] 

%%%  P(7)  =  PdotR,  rocket  roll  accel  (rad/s/s)  [f  (t)  ] 

%%%  Aero  Coefficients /Derivatives (functions  of  Mach  or  time): 

%%%  C(l)  =  CdO,zero  alpha  drag  coefficient 

%%%  C(2)  =  CNa, normal  force  /rad  alpha  (lift  curve  slope) 

%%%  C(3)  =  CMq,  pitch  damping  per  rad/sec  pitch  rate 

%%%  C(4)  =  Clp,  roll  damp  per  rad/sec  roll  rate  (HRGK) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  Related  Quantities 


%%% 

alpha. 

Angle  of  Attack 

%%% 

beta. 

Sideslip  Angle 

%%% 

Tz, 

Resultant  Vert,  Thrust  (up,  -z  dir.) 

%%% 

Ty, 

Reultant  Side  Thrust  (y  direction) 

%%% 

vt. 

Missile  Speed 

%%% 

Tertp, 

Air  Temperature  (K)  [f (altitude) ] 

%%% 

Pres, 

Air  Pressure  (Pa)  [f (altitude) ] 

%%% 

Rho, 

Air  Density  (Kg/m^3)  [f  (altitude)  ] 

%%% 

gbar. 

Dynamic  Pressure  =  1/2  rho(h)  Vt'^2 

%%% 

li. 

missile  length  (m) 

%%% 

d. 

missile  diameter  (m)  (ref  length) 

%%% 

S=pi*d''2/4, 

missile  cross-sectional  area (m^2) (ref  area) 

%%% 

mu. 

rocket  to  kit  coefficient  of  friction 

%%%  Aerodynamic  Forces 

%%%  CN  =  CNa*alpha  =>  N=CN*qbar*S+Tz  Normal  Force  («Lift) 

%%%  CY  =  -CNa*beta  =>  Y=CY*qbar*S+Ty  Sideslip  Force 

%%%  CD  =  CdO  +  CN^2/CNa  =>  D=CD*qbar*S  Drag  Force 

%%%  Aerodynamic  Moments 
%%%  CLbar  =  Clp*P*d/ (2*U) , 

%%%  =>  Lbar  =  CLbar* qbar*S*d  Kit  Roll  Moment 

%%%  CM  =  CNa*alpha(Xcg-Xcp)  +  CMq*Q*d/ (2*U) 

%%%  =>  M  =  CM*qbar*S*d  +  Tz(Xcg-Xt)  Pitch  Moment 

%%%  CNbar  =  CNa*beta(Xcg-Xcp)  +  CMq*Q*d/ (2*U) 

%%%  =>  Nbar  =  CNbar *qbar*S*d  +  Ty(Xcg-Xt)  Yaw  Moment 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  Declare  Global  Variables 

global  gO  TempO  PresO  B_atm  R_air  S  d  Sd  Xt  mu; 
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%%%  Define  State  Variables  from  Inputs 
X  =  U0(l:14) ; 

%%%  Body-“Axes  Velocities 

U  =  x(l)  ;  %  speed  along  missile  axis,  x 

V  =  x(2)  ;  %  speed  of  sideslip,  y  direction 

W  =  x(3);  %  speed  along  z  direction  ("Down") 

%%%  Body-Axes  Angular  Velocities 
P  =  x(4) ;  %  roll  rate 

Q  =  x(5) ;  %  pitch  rate 

R  =  x(6) ;  %  yaw  rate 

%%%  Euler  Angles 

phi  =  x(7);  %  roll 

theta  =  X ( 8 ) ;  %  pitch 

psi  =  x(9);  %  yaw 

%%%  Rocket  Rolling 

Pr  =  x(10);  %  rocket  roll  rate 

%%%  Location  Variables 

Pn  =  x(ll);  %  Position  North  of  (0,0,0) 

Pe  =  x{12);  %  Position  East  of  (0,0,0) 

h  =  X  ( 13 ) ;  %  Height 

%%%  Define  Control  Variables  from  Inputs 
u  =  uo(14:17) ; 

T1  =  u(l);  %  thrust  from  thruster  #1 

T2  =  u(2);  %  thrust  from  thruster  #2 

T3  =  u(3);  %  thrust  from  thruster  #3 

T  =  u(4);  %  thrust  through  the  missile  long,  axis 

%%%  Other  Missile  Parameter  Inputs  (time  or  Mach  dependent)  : 

Param  =  uo(18:24) ; 

mass  =  Param(l)  ;  %missile  mass  (kg) 

Xcg  =  Param (2);  %missile  Cg  (m  from  nose) 

Xcp  =  Param (3);  %itdssile  CP  (m  from  nose) 

JxR  =  Param (4);  %rocket  M  of  I  about  long  axis  (kg  m^2) 

JxK  =  Param(5);  %kit  M  of  I  about  long  axis  (kg  m^2) 

Jy  =  Param(6);  %missile  MofI  about  y  or  z  axes  (kg  m^2) 

PdR  =  Param(7);  %rocket  roll  acceleration  (rad/s/s) 

%%%  Define  Mach  Number  Dependent  Coefficients /Derivatives 
C  =  uo(25:28) ; 

CdO  =  C(l);  %zero  alpha  drag  coefficient 

CNa  =  C(2);  %normal  force  /rad  alpha  (lift  curve  slope) 

CMq  =  C(3);  %pitch  daitping  per  rad/sec  pitch  rate 

dp  =  C(4);  %roll  danping  (kit  only)  per  rad/sec  roll 
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%%%%%%%  Define  the  Forces  and  Moments 


%%%  Related  Quantities 
%  speed  quantities 
Vxz2  =  U^2  +  W^2; 

Vt2  =  Vxz2  +  V^2; 

Vt  =  sqrt(Vt2); 

%  atmospherics 

tertp  =  TenpO-B_atm*h; 

pres  =  PresO*(l-B__atm*h/TertpO)'"(gO/R_air/B_atm)  ; 
rho  =  pres/ {R_air* temp)  ; 
qbar  =  .5*rho*Vt2; 

%  control  thrusts 

Ty  =  Tl*cos(phi)  +  ... 

T2*cos  (phi+2/3*pi)  +  T3*cos  (phi+4/3*pi)  ; 

Tz  =  Tl*sin(phi)  +  ... 

T2*sin(phi+2/3*pi)  +  T3*sin(phi+4/3*pi)  ; 

%  angles  of  attack 
alpha  =  atan2(W,U); 
beta  =  atan2  (V,  sqrt  ( Vxz2 ) ) ; 


%%%  Forces 

CN  =  CNa*alpha; 

Fz  =  -‘qbar*S*CN  -  Tz; 
CY  =  -CNa*beta; 

Fy  =  qbar*S*CY  +  Ty; 
CD  =  CdO  + 

Fx  =  T  -  qbar*S*CD; 


%  normal  force  coefficient 
%  normal  force 
%  side  force  coefficient 
%  side  force 
drag  force  coefficient 
%  drag  force 


(CN^2+Cy^2)/CNa;% 


%%%  Moments 

Cl  =  Clp*P*d/ (2*U+ ,  001) ;  %  roll  .moment  coefficient 

Libar  =  qbar*Sd*Cl  +  mu*JxR*PdR;  %  roll  moment  on  kit 
CM  =  CNa*alpha*  (Xcg-Xcp)  +  ...  %pitch  moment  coefficient 
CMq*Q*d/(2*U+.001) ; 

M  =  qbar*Sd*CM  +  Tz*  (Xcg-Xt) ;  %  pitch  moment 

Cn  =  -CNa*beta*  (Xcg-Xcp)  +  ...  %  yaw  moment  coefficient 
CMq*R*d/(2*U+.001) ; 

Nbar  =  qbar*Sd*Cn  +  Ty*  (Xcg-Xt) ;  %  yaw  moment 


85 


%%%  Conpute  the  Time  Derivatives  from  Flat-Earth  Equations 
%%%  Force  Equations 

Ud  =  R*V  -  Q*W  -  gO*sin( theta)  +  Fx/mass; 

Vd  =-R*U  +  P*W  +  gO*sin (phi)  *cos  (theta)  +  Fy/mass; 

Wd  =  Q*U  -  P*V  +  gO*cos  (phi)  *cos  (theta)  +  Fz/mass; 

%%%  Moment  Equations 
Pd  =  Lbar/JxK; 

Prd  =  PdR; 

Qd  =  (l-JxK/Jy)*P*R  +  M/Jy  -  ( JxR/Jy)  *R*Pr; 

Rd  =  (JxK/Jy-l)*P*Q  +  Nbar/Jy  +  ( JxR/Jy)  *Q*Pr; 


%%%  Kinematic  Equations 

phid  =  P  +  tan(theta)  *  (Q*sin(phi)  +  R*cos(phi)); 

thetad  =  Q*cos(phi)  -  R*sin(phi); 

psid  =  (Q*sin(phi)  +R*cos  (phi) )  /cos (theta)  ; 

%%%  Navigation  Equations 

Pnd  =  U*cos  (theta)  *cos  (psi)  +  ... 

V*  (sin  (phi)  *sin(  theta)  *cos  (psi)  -  cos  (phi)  *s  in  (psi) )  +... 
W*  (cos  (phi)  *sin(  theta)  *cos  (psi)  +  sin(phi)  *sin(psi) )  ; 

Ped  =  U*cos (theta) *sin (psi)  +  ... 

V*  (sin  (phi)  *  sin  (theta )  *sin  (psi)  +  cos  (phi)  *cos  (psi) )  +... 
W*  (cos  (phi)  *sin{  theta)  *sin (psi)  -  sin(phi)  *cos  (psi) )  ; 
hd  =  U*sin(theta)  -  V*sin  (phi  )*  cos  (theta)  -  ... 

W*cos (phi) *cos (theta) ; 


%%%%%%%%%%  Define  the  output  vector 
%%%  Body-Axes  Accelerations 

xo{l)  =  Ud;  %  acc  along  missile  axis,x  direction 

xo{2)  =  Vd;  %  acc  of  sideslip,  y  direction 

xo{3)  =  Wd;  %  acc  along  z  direction  ("Down") 

%%%  Body-Axes  Angular  Accelerations 


xo(4)  =  Pd; 
xo(5)  =  Qd; 
xo(6)  =  Rd; 

%%%  Euler  Angle  Rates 
xo{7)  =  phid; 
xo(8)  =  thetad; 
xo(9)  =  psid; 

%%%  Rocket  Rolling  Equations 
xo(lO)  =  Prd; 

%%%  Location  Variable  Rates 
xo(ll)  =  Pnd; 
xo(12)  =  Ped; 
xo(13)  =  hd; 


%  roll  acc 
%  pitch  acc 
%  yaw  acc 

%  d  phi/dt 
%  d  theta /dt 
%  d  psi/dt 

%  rocket  spin  acceleration 

%  Position  North  of  (0,0,0)  rate 
%  Position  East  of  (0,0,0)  rate 
%  Height  rate 
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4.  Line  of  Sight 

fmiction  yo  =  Lineof Sight  (uo) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  This  Function  computes  the  seeker  frame  LOS, 

%%%  from  the  tangent  plane  coordinates 
%%%  W.  Mark  Wonnacott,  May  1997 
%%%  Naval  Postgraduate  School 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%% 

INPUTS 

%%% 

Position  Differences  Xt-Xm 

%%% 

u(l)  =  N, 

delta  north 

%%% 

u(2)  =  E, 

delta  east 

%%% 

u(3)  =  H, 

delta  height 

%%% 

Euler  Angles 

%%% 

u(4)  =  phi, 

seeker  roll  relative  to  vertical  down 

%%% 

u(5)  =  theta, 

rocket  pitch  relative  to  horizontal 

%%% 

u(6)  =  psi, 

rocket  yaw  relative  to  north 

%%% 

OUTPUTS  Seeker  LOS 

%%% 

y(l)  =  alpha. 

off-boresight  angle 

%%% 

y(2)  =  beta. 

angle  from  seeker's  ref.  roll  position  (down) 

%%% 

y(3)  =  R,  slant  range  to  target 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  Define  Tangent  Plane  LOS  Offsets 

N  =  uo(l);  %  north  distance  to  target 

E  =  uo(2);  %  east  distance  to  target 

D  =  uo(3);  %  down  distance  to  target 

%%%  Define  Rocket  Euler  Angles 

phi  =  uo(4);  %  seeker  roll  relative  to  vertical  down 

theta  =  uo(5);  %  rocket  pitch  relative  to  horizontal 

psi  =  uo(6)  ;  %  rocket  yaw  relative  to  north 

%%%  Transformation  Matrices  (NED  to  XYZ-body) 

Cl  =  [cos  (psi)  sin(psi)  0; 

-sin (psi)  cos  (psi)  0; 

0  0  1]  ; 

C2  =  [cos (theta)  0  -s in ( theta ) ; 

0  1  0; 

sin ( theta)  0  cos ( theta) ] ; 

C3  =  [  1  0  0; 

0  cos  (phi)  sin  (phi )  ; 

0  -sin (phi)  cos(phi)]; 

X=C3*C2*C1*[N;E;D]; 

Offbore  =  atan2  (sqrt  (X(2)  ^2+X(3 )  ^2)  ,X(1) )  ; 
pointing=  atan2  (X  ( 3 ) ,  X  ( 2 ) )  ; 

R  =  sqrt(X(l)^2+X(2)^2+X(3)'"2); 

yo(l)  =  Offbore;  yo(2)  =  pointing;  yo(3)  =  R; 
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5.  Mach  Number 

fxinction  M  =  Mach(uo) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  This  fimction  catputes  the  Mach  NiJinber  as  a  function  of 
%%%  conponent  velocities  and  altitude 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  W.  Mark  Wonnacott,  M^ch  1997 
%%%  Naval  Postgraduate  School 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%% 

INPUTS: 

%%% 

Body-Axes  Velocities 

%%% 

x(l)  =  U, 

speed  along  missile  axis,  x  direction 

%%% 

x(2)  =  V, 

speed  of  sideslip,  y  direction 

%%% 

x(3)  =  W, 

speed  along  z  direction  { "Down" ) 

%%% 

x(13)=  h, 

altitude  (m) 

%%% 

Related  Quantities 

%%% 

Vt, 

Missile  Speed  (m/s) 

%%% 

Terrp, 

Air  Teitperature  (K)  as  function  of  altitude 

%%% 

Pres, 

Air  Pressure  (Pa)  as  fimction  of  altitude 

%%% 

Rho, 

Air  Density  (Kg/m^3)  as  function  of  altitude 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  Declare  Global  Variables 

global  gO  TenpO  PresO  B__atm  R_air  S  d  Sd  Xt; 


%%%  Define  Variables  from  Inputs 
%%%  Bocfy-Axes  Velocities 

U  =  uo(l);  %  speed  along  missile  axis,  x  direction 

V  =  uo{2);  %  speed  of  sideslip,  y  direction 

W  =  uo(3);  %  speed  along  z  direction  ("Down") 

h  =  uo  ( 13 ) ;  %  height 


Vt  =  sqrt(U^2  +  V"^2  +  W^2); 


%  atmospherics 
k  =  1.4; 

tenp  =  TerrpO-B_atm*h; 

pres  =  PresO*  (1-B_atm*h/Teirp0)  ^  {gO/R_air/B_atm)  ; 
rho  =  pres/ (R_air*terrp) ; 


c  =  sqrt  (k*R_air*tertp)  ; 
M  =  Vt/c; 
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